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DIFFUSION APPROXIMATION FOR AN OVERLOADED 
X MODEL VIA AN AVERAGING PRINCIPLE 

By Ohad Perry* and Ward Whitt^ 
Northwestern University and Columbia University 

In previous papers we developed a deterministic fluid approx- 
imation for an overloaded Markovian queueing system having two 
customer classes and two service pools, known in the call-center lit- 
erature as the X model. The system uses the fixed-queue-ratio-with- 
thresholds (FQR-T) control, which we proposed in a recent paper 
as a way for one service system to help another in face of an unex- 
pected overload. Under FQR-T, customers are served by their own 
service pool until a threshold is exceeded. Then, one-way sharing is 
activated with customers from one class allowed to be served in both 
pools. The control aims to keep the two queues at a pre-specified fixed 
ratio. We supported the fluid approximation by establishing a many- 
server heavy-traffic functional weak law of large numbers (FWLLN) 
involving an averaging principle. In this paper we develop a refined 
diffusion approximation for the same model based on a many-server 
heavy-traffic functional central limit theorem (FCLT). 



1. Introduction. In this paper we establish a many-server heavy-traffic 
functional central limit theorem (FCLT) for an overloaded large-scale Marko- 
vian queueing system having two classes and two service pools, known as the 
X model [7], using the fixed- queue-ratio with thresholds (FQR-T) routing, 
which we proposed in [21]. 

In particular, we consider a system in which each class has its own des- 
ignated service pool, but with all agents, in both pools, capable of serving 
customers from both classes. The control aims to prevent sharing of cus- 
tomers (i.e., sending customers from one class to be served at the other class 
pool) when both classes are normally loaded, and activate sharing when the 
system unexpectedly experiences an overloaded, due to an unforseen shift 
in the arrival rates. 

When sharing is taking place, the control aims at keeping a pre-specified 
fixed ratio between the two queues. This ratio is chosen according to a 
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deterministic ("fluid") optimization problem; see §5.3 in [21], where it is also 
shown that sharing should not be allowed at both directions simultaneously, 
i.e., at any time there should be at most one pool working with both classes. 
In general there are two different ratios: If class 1 is overloaded, then an 
optimal ratio r\2 should hold between the queues. If class 2 is overloaded, 
then an optimal ratio r2,i should hold between the queues. In [23] we showed 
that the FQR-T control achieves the objectives above, asymptotically in the 
fluid limit. Moreover, the FQR-T control produces a tractable fluid limit. 
Here we establish a refined stochastic limit. 

The FQR-T control here is a modification of the FQR control (without 
the thresholds), which is a special case of the queue- and-idleness ratio (QIR) 
controls suggested in [11]. These QIR and FQR controls were analyzed in [10, 
11, 12] for critically loaded systems, operating in the quality and efficiency 
driven (QED) many-server heavy-traffic regime; see [8, 13]. Heavy-traffic 
limits for networks having cyclic graphs, such as the X model, were obtained 
under the condition that the service rates are class or pool dependent; see 
Theorem 3.1 in [11]. In general, when the service rate depends on both the 
class and the pool, FQR can perform badly in cyclic networks, creating 
severe congestion even if each pool is not congested by itself; see §4.1 in [21] 
and §EC2 in [22]. 

We suggested the FQR-T control in [21], and analyzed the X model using a 
heuristic stationary fluid approximation. In [22] we determined the transient 
behavior of that same fluid model, based on a stochastic averaging principle 
(AP), but that AP was introduced there as a heuristic engineering principle, 
supported only by simulation. The purpose of our subsequent papers [23, 24] 
was to establish key mathematical properties of the fluid model, expressed 
as an ordinary differential equation (ODE), and show that the fluid model, 
heuristically derived in [21, 22], arises as the many-server heavy traffic limit 
of a sequence of X models in the many-server efficiency driven (ED) regime. 
That FWLLN is challenging, because the fluid limit depends critically on 
the AP. For each n, the system evolves as a 6-dimensional continuous-time 
Markov chain (CTMC), but there is (a somewhat complicated) statistical 
regularity associated with the many-server heavy-traffic limit. In particular, 
the limiting fluid approximation is a deterministic function characterized 
by an ODE (and an initial condition), which is driven by the time- varying 
instantaneous average behavior of a family of fast-time-scale stochastic pro- 
cesses (FTSP's), which produces the AP. See §1.3 of [24] for a discussion of 
the literature on AP's; notable contributions in the queueing literature are 
[4, 15]. See [6] for a FCLT involving an AP, building on [15]. 

We now build on the FWLLN and the AP to describe the distribution of 
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the stochastic fluctuations about the fluid path; i.e., we establish the cor- 
responding FCLT, which is Theorem 4.1 here. There is technical novelty in 
properly treating the FTSP's alluded to above. The limit process involves 
an independent Brownian motion term with deterministic time scaling in- 
volving the asymptotic variance of the FTSP; see §4.1 and L2, I, 72 and 73 
in Theorem 4.1. A key step in establishing the main result - the FCLT in 
Theorem 4.1 - is a FCLT for the family of FTSP's, Theorem 6.1, which is 
of independent interest. This challenging step proves a FCLT for a sequence 
of CTMC's having time-varying parameters depending on the fluid limit. 

From an engineering perspective, Corollary 4.1 is especially useful for un- 
derstanding the performance of the FQR-T control. It describes the stochastic- 
process limit once the fluid has stabilized (i.e. when the fluid is stationary). 
With a constant fluid state, the key limit process becomes the well-studied 
bivariate Ornstein-Uhlenbeck (BOU) process, which has a Gaussian distri- 
bution for each t; see Corollary 4.1 below. Consequently, the approximating 
steady-state distribution during the overload is a Gaussian distribution, with 
mean values equal to the stationary fluid point in Theorem 3.2 multiplied 
by n, and variance and covariance terms in (4.16) multiplied by y/n. 

The FCLT extension is essential for truly understanding the system per- 
formance under overloads, because the actual performance is not nearly 
deterministic, as described by the fluid approximation, unless the scale is 
extremely large. This phenomenon is well illustrated by the example here 
in §11. For that example, the standard deviations of the queue lengths are 
about equal to (half of) the mean queue lengths when the number of servers 
in each pool is 25 (100). 

Here is how the paper is organized: After preliminaries in §2, we briefly 
state the FWLLN and the associated WLLN for the stationary distribu- 
tions in §3. We state the FCLT and our other main results in §4. We prove 
the FCLT in §5 except for Lemma 5.6, establishing joint convergence of the 
driving processes. We give the proof of Lemma 5.6 in §6 except for two key 
supporting results. One is an approximation of the queue-difference pro- 
cess by associated frozen processes (having states frozen at initial values) 
in Lemma 6.2 and the other is a FCLT for the FTSP with time-varying 
parameter state function in Theorem 6.1. We prove Lemma 6.2 in §7 and 
Theorem 6.1 in §8. Our proof of Lemma 8.4 to prove Theorem 6.1 exploits 
the martingale FCLT for triangular arrays. We state these supporting mar- 
tingale results in §9. We then prove all remaining lemmas in §10. Finally, 
we evaluate the quality of the approximations by making comparisons with 
simulations in §11. 
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2. Preliminaries. 

2.1. Notation. Let R, Z and N denote the real numbers, integers and 
nonnegative integers, respectively. Let = denote equality by definition. For 
a subinterval / of [0,oo), let V = V(I) = V(I, R) be the space of all right- 
continuous R-valued functions on / with limits from the left everywhere, 
endowed with the familiar Skorohod J\ topology [32]. Let C be the subset 
of continuous functions in T>. Let a subscript k appended to one of these 
spaces denote the set of all ^-dimensional vectors with components from the 
space, endowed with the corresponding product topology, e.g., R& and T>^. 

Let dj 1 denote a metric on T>k(I) inducing the convergence. Since we 
will be considering continuous limits, the topology is equivalent to uniform 
convergence on compact subintervals of /. Let e be the identity function 
in D = T>i; i.e., e(t) = t, t & I. Let o be the composition function, i.e., 
(.x o y)(t) = x(y(t)). Let denote convergence in distribution [32]. 

We use the familiar big-O and small-o notation for deterministic functions: 
For two real functions / and g, we write 

f(x) = 0(g(x)) whenever limsup \f{x)/g(x)\ < oo, 

X— >OD 

f{x) = o{g(x)) whenever limsup \f{x)/g{x)\ = 0. 

x— >oo 

For a function x : [0, oo) — > R and < t < oo, let ||x||t = sup 0<s<i 

For a stochastic process Y = {Y(t) : t > 0} and a deterministic function 
/: [0, oo) ->• [0, oo), we say that Y is op(f(t)) if ||Y|| t //(t) as t ->• oo. 

For a sequence of stochastic processes or random variables, {Y n : n > 1}, 
we denote its fluid-scaled version by Y n = Y n /n. We let Y n = Y n /^/n be 
the -y/n-scaled processes without the centering about the fluid limit, and Y n 
denote the diffusion-scaled processes centered about the fluid limit, as in 
(4.5) below. 

2.2. A Sequence of Overloaded Markovian X Models. We consider a se- 
quence of overloaded Markovian X models, indexed by superscript n. There 
are two customer classes and two service pools. We are looking at these 
models during the overload incident, after the arrival rates have changed. 
The arrival rates are considered fixed, but the system is typically not yet in 
its new steady-state during the overload (assuming that the overload would 
persist). For each n and i = 1, 2, there is a class-i Poisson arrival process with 
rate A" . Customers have limited patience, and may abandon when waiting 
in queue. The times to abandon are i.i.d. exponential variables with rate 
9i for each class-i customer in queue. Service pool j has ml- homogeneous 
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agents (servers). Service times of class- i customers by pool-j agents are mu- 
tually independent and exponentially distributed with rate mj, i,j = 1,2. 
The abandonment and service rates are independent of n. 

Since we are considering an overload incident, we will scale to achieve an 
efficiency-driven (ED) many-server heavy-traffic regime. 

Assumption 1. (many-server heavy-traffic scaling) For Aj,m,, > 0, i = 
1,2, 

— - — = > ana -= > U as n — > oo. 

We could instead obtain a modified, more general, FCLT if there were 
nondegenerate limits in Assumption 1, but we consider our choice natural, 
because the system operates in an overload regime. (For the FWLLN, it is 
sufficient that Xf/n — > A.; and mf/n — > mi as n — > oo, i = 1, 2.) 

Let 

= and tfES ^-^™^ , < = 1,2, 

where, for y S R, y + = max{0,y}. Then pi is the traffic intensity for pool i 
and qf is the stationary class-i fluid-limit queue, when both pools operate 
independently. We say that pool i is overloaded if pi > 1. However, with 
sharing allowed, pool i can be overloaded even if pi < 1 provided that 
enough class j customers are routed to be served there, j ^ i. The next 
assumption makes precise our notion of system overload. 

Assumption 2. (system overload, with class 1 more overloaded) 
The rates in the system are such that (1) 9iqf > ^i,2 m 2(l — P2) + and (2) 

Clearly, p\ > 1 by Condition (1), so that class 1 is overloaded. However, 
Condition (1) also ensures that pool 2 is overloaded if sharing is taking place. 
That is so because, even if p 2 < 1, there is not enough extra service capacity 
in pool 2 to take care of all the class-1 customers that pool 1 cannot serve. 
Condition (2) in the assumption implies that even if pool 2 is overloaded by 
itself (i.e., if p 2 > 1), then class 1 is the one that should receive help from 
pool 2. 

2.3. The FQR-T Control. We now describe the FQR-T control for each 
system n. The purpose of the FQR-T control is: (i) to prevent sharing under 
normal loads, (ii) to activate sharing as soon as an overload incident begins, 
and (iii) to keep close to the desired ratio between the two queues, making 
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sure that sharing takes place in the needed direction only. The control is 
based on two positive thresholds, k™ 2 an d fc^i) an d the two ratio parameters 
discussed above, r\ >2 arid r 2i \, which satisfy r± i2 > r 2,i; see Proposition EC. 2 
and Equation (ECU) in [21]. 

Let Qf(i) be the number of customers in the class-i queue and let ZfAt) 
be the number of class-i customers being served in service pool j, at time t, 
i,j = 1,2 (in the n th system). The FQR-T routing is based on the queue- 
difference stochastic processes 

D't 2 (t)^Q r Kt)-k'i 2 -r lt2 Q^t), and 
DIS) = r 2 ,iQ2(t) - k$ tl - Q?(t), t > 0. 

As long as Z?" 2 (^) — an d DQiit) < 0, no sharing of customers is allowed, 
i.e., a server in pool j takes only class j customers, j = 1, 2. It follows from 
[8] that thresholds of order larger than 0(yfn) will prevent sharing in such 
circumstances, asymptotically as n — > oo. Once one of the queue-difference 
processes in (2.1) becomes strictly positive (so that one of the thresholds is 
crossed) sharing is initiated. It follows from the Corollary 2.1 in [33], that 
thresholds of size o(n) will detect an overload relatively quickly (instantly, 
asymptotically as n — > oo). We thus choose the thresholds according to the 
following assumption. 

Assumption 3. (scaling of the thresholds) For ki 2 ,k 2t \ > and a se- 
quence of positive numbers {c n : n > 1}, where c n /n — > and c n /y/n — > oo 
as n — > oo, 

fe" 2 /cn — > &i,2 and k,2i/c Tl — > k 2> \ as n — > oo. 

Finally, only one-way sharing is allowed at any time. For example, a newly 
available pool-2 agent at time t serves a class-1 customer if D™ 2 {t) > 0, 
provided no class-2 customers are served in pool 1 at that same time t; 
otherwise he serves a class-2 customer. 

2.4. Dimension Reduction. For the X model operating under FQR-T, 
the six-dimensional process 

(r, r)\ -y-n / r\Tl /~\rl ryn ryn ryn ryn \ 

l z - z J -^-6 = Wl i ^62i ^l,2i -^2,1) ^2,2) 

is a CTMC for each n > 1. However, there is an important dimension re- 
duction established in §6 of [24]. It was shown, under the assumptions above 
and with appropriate initial conditions, that asymptotically the two service 
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pools remain fully occupied with no pool-1 servers serving class 2; i.e., for 
each T > 0, 

P(Z[ l 1 {t) = m'?,Zl 1 {t) = 0,Z r 1 \ 2 + Zl 2 = m%,0<t<T)^l as n -»• oo. 
Thus, the system is characterized by an essentially three-dimensional process 

(2.3) X%>* ee (Q?,Q^m?,^ 2) 0,m^ - Z^), 
having the vector of essential components 

(2.4) X"ee(Q?,Q« Z£ 2 ), 

whose evolution is directly specified, and will be specified here in Theorem 
5.1. Theorem 5.1 concludes that Xq '* and X£ are asymptotically equivalent, 
so that X n is sufficient to characterize the FWLLN and, in turn, to prove 
the FCLT. That implies that X% => x 6 in £> 6 if and only if X n => x in T> 3 as 
n — > oo, with x(t) E S = [0, oo) 2 x [0, 771,2], for t > 0; see Theorem 3.1 below. 
We thus restrict attention to the space S. 

2.5. The Fast- Time- Scale Process. Given that the system is overloaded 
with class 1 needing help from pool 2, as determined by Assumptions 1 
and 2, the FQR-T control is driven by the process 2 i n (2-1)- Since the 
queue lengths are asymptotically of order 0(n), the queue-difference process 
D™ 2 has transitions at rate 0(n). However, Theorem 4.4 in [24] shows that, 
under regularity conditions, the sequence {D™ 2 (t) '■ n > 1} is stochastically 
bounded in M, so that the difference process should be analyzed without any 
spatial scaling. On the other hand, Theorem 4.4 in [22] also shows that this 
sequence is not P-tight. Thus, these difference processes do not converge to 
nondegenerate limits in T> as n — > 00 without spatial scaling. Nevertheless, 
both the FWLLN and FCLT depend heavily on the asymptotic behavior of 
functionals of that driving queue-difference process and on the analysis of a 
related family of fast time scale process (FTSP's). 

Fix t > and consider the time expanded queue- difference process 

(2.5) {D n e (T n , s) : s > 0} = {Z^ 2 (i + s/n) : s > 0}, 

where r n is a deterministic vector in R3, and we condition on X n (to) = 
T n . Theorem 4.3 in [24] shows, under the assumptions of the FWLLN in 
Theorem 3.1 below, that 

(2.6) {D T e l (T n , s) : s > 0} =► {£(7, s) : s > 0} in V as n -> 00 
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if r n — 7- 7 G § and D™(T n , 0) => D(j, 0) in R as n ->■ oo. The limit process 
-D(7, •) is the FTSP, an irreducible pure-jump (time homogeneous) Markov 
process having transition rates that are the limit of the instantaneous rates 
of £>i j2 (*o) at time t (given the state of the CTMC Xg(t )) } divided by n. 
Since the distribution of the FTSP is determined by 7, we obtain a different 
FTSP D(7, •) for each 7 G S, and thus for each t > 0. The name "FTSP" 
becomes clear when observing that it arises as the limit in (2.5) achieved by 
"slowing" time in the neighborhood of each time point to in D± 2 {to)- 

In this paper we will be concerned with the region in which the FTSP's 
are positive recurrent. In particular, we let A denote the set of all 7 G §, 
such that Z?(7, •) is positive recurrent, with .0(7,00) denoting a random 
variable distributed as the stationary distribution of the FTSP D(j, •). For 
each 7 G §, let 

(2.7) 7n, 2 ( 7 ) = P(£>(7,oo) >0). 

By Lemma 3.1 in [24], 717,2(7) is well defined for all 7 G S, but 0(7,-) is 
positive recurrent if and only if < 717,2(7) < 1. By Theorem 6.1 of [23], 

(2.8) A = {7 : < 711,2(7) < 1} = {7 : 5 + (j) < and £-(7) > 0}, 

where <5+ (7) and #-(7), respectively, are the constant drift rates in the pos- 
itive region {s : D(7, s) > 0} and the non-positive region {s : D(-f, s) < 0}. 

Both the FWLLN and the FCLT depend critically on distributional and 
topological characteristics of the FTSP's. A simplification is achieved by 
representing the FTSP as a quasi-birth- and- death (QBD) process, which 
can be done by assuming that 77,2 is rational. We refer to §6.2 in [23] for 
more details on the QBD representation of the FTSP, and to [18] for the 
general theory of QBD processes. See also Theorem 6.1 and Equation (7.2) 
in [23] for how the QBD representation simplifies the characterization of A, 
as well as §11 in [23], where an efficient algorithm for computing the fluid 
limit numerically is developed, based on that QBD representation. We thus 
make the following assumption. 

Assumption 4. {queue ratios parameters) The queue ratios 7*1,2 andr2 ) i 
are positive rational numbers. 

Since we are considering the case when sharing is taking place with class-1 
customers receiving help, we essentially need only consider 77 2, which we 
henceforth denote by r, i.e., r = ri,2- 
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3. The Fluid Limit. We now review the FWLLN for the process X§ 
in (2.2) and the WLLN for the associated sequence of stationary random 
variables (oo), established in [24]. For these, we assume that the fluid 
x(t) is in the set A, where the FTSP is positive recurrent. We conclude by 
reviewing a result stating that the fluid model eventually remains in A. 

3.1. The FWLLN. We now describe the fluid limit, i.e., the limit of X% 
for Xq in (2.2). The FWLLN requires an assumption about the initial con- 
ditions. In [24] we considered a (more general) version of the following. 

Assumption 5. Assume that 

P(Z£ 1 (0) = 0,Q?(0) > a n ,i = 1,2) = 1 for all n > 1, 
X n (0) =4> x(0) € A and D™ 2 (0) => L as n -too, 

where L is a finite random variable, x(0) is deterministic and {a n : n > 1} 
is a sequence of numbers satisfying a n /c n — >• a, < a < oo, for c n in 
Assumption 3. 

We note that in [24] x(0) was not necessarily in A. The following theorem 
is a version of the main result - Theorem 4.1 - in [24], adapted to our needs 
here. 

Theorem 3.1. (FWLLN) Under Assumptions 1-5, 

Xq xq in Vq([0,oo)) as n — > oo, 

for Xq in (2.2), where x§ = (%, Zij;i, j = 1,2), is a deterministic element 
of Cq, with zi t i = mi, z-i,\ = Oe and Z2,2 = fU2 — and x = (91,92,-^1,2) 
being the unique solution to the three-dimensional ODE 

qi(t) = Ai - mi/xi,! - 7ri j2 (a;(i)) [21,2(^^1,2 + ^2,2(^)^2,2] - ^i9i(*) 
(3.2) q 2 (t) = A 2 - (1 - 7ri )2 (x(t))) [^2,2(^)^2,2 + ^,2(^^1,2] - #292^) 
zi,2(t) = n 1>2 (x(t))z2,2(t)n2 l 2 ~ (1 - vri i2 (x(t)))zi i2 (t)^i, 2 , 

for Tt\2{x{t)) = P(D(x(t), 00) > 0) in (2.7). Moreover, there exists 5, < 
5 < 00, such that x(t) G A, so that < 7Ti j 2(x(t)) < 1 and qi(t) = rq2{t), 
for allte [0,5). 

Just as the routing of customers at each time t > in the prelimit is 
determined by whether D™ 2 (^) > or < 0, so also the instantaneous future 
evolution of the fluid limit x(t) at time t > 0, is determined by whether the 
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FTSP corresponding to x(t), D(x(t), •), is positive or nonpositive. However, 
that evolution is determined by the long-run average behavior of the FTSP 
corresponding to time t, i.e., by TTi^{x(t)), giving rise to the term "averaging 
principle". Loosely speaking, D" 2 (t) achieves a local steady state (the steady 
state of the FTSP) instantaneously as n — > oo, at each time t > 0. 

Observe that Theorem 3.1 concludes that, if x(t) G A, then q\(t) = rq2(t). 
(That part of the theorem follows from Theorem 4.4 in [24].) It follows that 
A C § is a two-dimensional subspace of S, so that we have SSC in the 
sense that the original six-dimensional process is a deterministic function of 
a two-dimensional process. More importantly for the FCLT, we will show 
that Q1(t) - k? >2 - rQ%(t) = o(y/n) for t G (ii,i 2 ) if x(t) G A over [ii,i 2 ), 
so the SSC to two dimensions holds in diffusion scale as well; see Theorem 
5.2 below. 

3.2. The Stationary Fluid Limit. Our main theorem here will be estab- 
lishing the FCLT about the fluid trajectory, given that the trajectory is in 
A. An important consequence will be the BOU limit when the fluid limit is 
stationary. Since the fluid limit of X n in (2.4) is the unique solution to the 
ODE (3.2), there is an immediate equivalence between stationarity of the 
fluid limit and stationarity of the dynamical system in (3.2), and we do not 
distinguish between the two. 

Definition 3.1. (fluid stationarity) A point x* G § is a stationary point 
of the unique solution x = {x(t) : t > 0} to the ODE (3.2) if x(0) = x* 
implies x = x*e. If x = x*e, then x is said to be stationary. 

Since the ODE is autonomous (i.e., time invariant), we can replace time 
with any t > in the definition 3.1. That is, if x(T) = x* for some T > 0, 
then x(t) = x* for all t > T. Time invariance also implies that x(t) is 
stationary at time t (x(t) = x*) if and only if x(t) = (qi(t), q2{t), ii,2 (*)) = 
(0,0,0); see §8 of [23]. 

There are several issues regarding stationarity, which we addressed in 
[23]. In advance, neither existence of a stationary point to the fluid limit nor 
uniqueness are immediate. Even if there exists a unique stationary point, 
it needs to be identified. Moreover, it must be shown that the fluid limit 
converges to a stationary point as t — > oo. (There are still other issues 
regarding stability of the dynamical system in (3.2), and we refer to §8.3 in 
[23] for a discussion.) Finally, the fluid limit of Xq in (2.2) is characterized by 
the fluid limit of the three-dimensional X n in (2.4), but that does not directly 
imply any relation between the stationary fluid limit and the stationary 
stochastic prelimit. 
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We now present the most relevant results for the FCLT regarding fluid 
stationarity. 

Theorem 3.2. (fluid stationarity) Under Assumptions 1-5, the following 
hold: 

(i) For each n, X^(t) =>■ Xg(oo) in R as t — > oo, with (oo) being 
the unique stationary distribution of the CTMC, and Xq (oo) =^ Xg in R as 
n — ?• oo /or 

(3.3) x* e = (q*,q2,rni,z* )2 ,0,m 2 - z{ i2 ), 
w/iere 

„ 6» 2 (Ai - mi/Ji i) - r0i(\ 2 - m 2 ^ 2 ,2) . 
z i,2 = a r~a A m 2> 

Al - mi/il,l - /il 2^1 2 , * A 2 - /U 2 ,2("l2 - 4 2 ) 

91 = ^ ~ an 92 = #2" ~ ■ 

(ii) x* = (ql,q 2 ,z^ 2 ) ^ s the unique stationary point of x, the unique so- 
lution to the ODE (3.2). 

(Hi) 711,2 (a;*) = P(D(x* ,00) > 0) = nl 2 , where D(x*,oo) is a random 
variable with the stationary distribution of the FTSP D(x*, •) and 

(3.4) 7T 12 



Ml.2^,2 + ( m 2 - 4,2)^2,2 ' 

(iv) x(t) — > x* as t — > 00 exponentially fast. 

Proof. Parts (i), (ii) and (Hi), and (iv), respectively, are covered by 
Theorem 4.2 in [24], §8 of [23] and Theorem 9.2 in [23]. Explicit exponential 
bounds on the rate of convergence to stationarity in (iv) are given in [23]. 
We now elaborate on (ii) and (Hi). First, if x* ^ A, then the fact that x* is a 
stationary point of x follows immediately from the fact that 7Ti j2 (x*) = or 
= 1. In that case, it is also easy to see that ir* 2 in (3.4) is equal to tt\ )2 (x*); 
see Corollary 8.1 in [23]. It is the unique stationary point by Theorem 8.1 
in [23]. The more challenging case, in which x* £ A and the existence of a 
stationary point is nontrivial, is proved in Theorem 8.2 in [23]. □ 

3.3. Eventually Remaining in the Set where the FTSP is Positive Recur- 
rent. The FCLT will be stated under the assumption that the associated 
fluid limit lies in the set A. Thus we now explain why this makes sense and 
introduce an additional assumption. 
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Note that x* G in (3.3) is completely characterized by x*, which involves only 
the rates in the system, and does not require any knowledge of the transient 
fluid limit or the initial condition. (In particular, SSC to three dimensions 
holds for the WLLN of the stationary distributions.) Simple algebra shows 
that if < z\ 2 < m 2 , then q\ = rq%. Together with (2.8) and (3.4) we see 
that x* G A if and only if < z\ 2 < ?7i2 . It follows from Assumption 2 and 
(3.3) that z\ 2 > (see also Corollary 8.2 in [23]), so that, under Assumption 
2, 

(3.5) x* G A if and only if z* 2 < ni2- 

The next theorem, which follows from Theorem 10.2 in [23], shows that there 
is not much loss in assuming that the limit x lies entirely in A whenever 
x* G A. 

Theorem 3.3. If x* G A then there exists Ta < oo such that x(t) G A 
for all t > Ta- 

Since we are interested in the case x* G A, which is the main case, as is 
clear from (3.5), we make the following assumption 

Assumption 6. For all t>0, x(t) g A. 

Assumption 6 is not essential for our results; we make it only for simplicity 
of the exposition. Without this assumption, the FCLT can be proved over 
a finite interval over which a; £ A. In applications, the fluid limit is likely 
to hit A immediately after the overload begins, and remain in A thereafter; 
see §11.3 in [23]. 

4. The Main Results. In preparation for the FCLT, we indicate how 
the limit is affected by the FTSP in §4.1. We then state the main FCLT and 
important corollaries in §4.2 and §4.3. We conclude in §4.4 by indicating 
how the results simplify in the special case r = n j = 1, where FQR reduces 
to serving the longer queue. 

4.1. The Role of the FTSP 's in the Stochastic Limit. Just as the limiting 
ODE in (3.2) arising in the FWLLN depends on the FTSP's D(j, •) (through 
the probability TTi t 2(x(t))), so too the stochastic limit process arising in the 
FCLT refinement depends on these same FTSP's. Since the FTSP -D(7, •) 
depending on the state 7 is a positive recurrent QBD under the assumption 
that 7 G A, the stochastic refinement depends on the asymptotic variability 
of the FTSP. In particular, since the FTSP D(j, •) is a regenerative process 
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(which can be represented as a QBD whenever the ratio r is rational), the as- 
sociated cumulative process obtained by integrating the indicator functions 
1 {D( 1 ,.s)>o} obeys a FCLT; i.e., 

rnt 

(4.1) n- 1 ' 2 (l{D( 7l -)>o}-7n.2(7)) ds^B(a\ 7 )t) in V 

Jo 

as n — > oo, where B is a standard Brownian motion (BM) for each 7 € A. 

The constant cr 2 (7) appearing inside the BM on the right in (4.1) is often 
called the asymptotic variance (see [3, 9, 31]) of the regenerative process 
D(j, s) (and the function / with f(D{^f, s)) = 1{_d( 7jS )>o})- For each 76 A, 
it is defined as the limit 

(4.2) a 2 (j) = fim jVar (J l {D{l)S>0} ds^\ . 

In our QBD context, the underlying regenerative cycles can be determined 
by successive visits of D(j, •) to any fixed state. The asymptotic behavior 
is determined by the random length of a cycle, and either the random 

integral over a cycle, ^(7), or the random centered integral over a cycle, 
^(7), where 
(4.3) 

Y(l)= 1{D( 7 , S )>0} d« and Y(-y) = / (1{d( 7 ,s)>o} - ^1,2(7)) ds. 



The key asymptotic quantities here can be expressed in terms of the means 
of the first two variables and the variance of Y{^) via 

4.4 7Ti, 2 7 = F 7 \i and a 2 7 = v \ ' 

see [3, 9]. Of course, ^(7) = ^(7) — ^i^il^il), so that Var(Y(^/)) can be 
expressed in terms the means, variances and the covariance of the variables 
t(7) and ^(7), where < ^(7) < t(7) w.p.l. Here we have strong regularity, 
with the random variable ^(7) having a finite moment generating function 
and all these quantities being continuous functions of the state 7, by virtue 
of Lemma A. 5 of [24]. 

4.2. The FCLT. Let Af(t) count the number of class-i customer arrivals, 
let Sfj(t) count the number of service completions of class- i customers by 
agents in pool j, an let U™{t) count the number of class-i customers to 
abandon from queue, all in model n during the time interval [0, t]. Let D™ 2 (i) 
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be the queue-difference process in (2.1) and let Q™(t) = Qi(t) + Q^it), all 
at time t. Let p\ = r/(l + r) and p 2 = 1 — p\ = 1/(1 + r), where r = T\i- 
For t > and i, j = 1, 2, let the diffusion-scaled processes be 





- n\i{t) 




fn 




- nzij(t) 




\/n 


Q?(t) ■ 


- nqi(t) 




fn 


Q?(t) ■ 


- nq s (t) 


Jn 



U?(t)-n0ifiqi(8)d8 



n 



n 



(4.5) ^ )e m^), gn (t) s gm 



^i n 2 (*) 



72 

/ n (i) ee \Zti / (1 {D » 2 ( S ) >0 } - ir h2 {x(s)))ds, t > 0, 
J 

where a; = (gi, q2, z\ 2) is the customary three-dimensional representation 
of the fluid limit, 2^1 = m\e, z 2 ,\ = 0e 22,2 = i^e — zip, q s = q\ + 
^2 and ^1,2(^(5)) ee P(D(x(s), 00) > 0), with D(a;(s),oo) being a random 
variable with the steady-state distribution of the FTSP {D(x(s),t) : t > 0} 
associated with the fluid limit at time s. 

Here is the main result of this paper: the FCLT for the overloaded X 
model operating under FQR-T. Since the limit is clearly a Markov process 
with continuous sample paths, it is by definition a diffusion process. Most 
of the rest of the paper is devoted to its proof. 

Theorem 4.1. (FCLT) //, in addition to Assumptions 1-6, 

(4.6) ((^(0), ^ 2 (0)) (4(0), £1,2(0)) G R 2 as n -> 00, 
i/ien, /or i, j = 1, 2, 

(4.7) (A?, U?, S^,D n , I n , Ql, Q™, =► (ii, ft, Sij,D, I, Q it Q„ 4 

in P17, where the processes depending on n on the left are defined in (4.5) 
and the limit process has continuous paths w.p.l. The initial 10- dimensional 
component (Ai,Ui, Sij, D, I) is a vector of independent Brownian motions, 
time scaled by increasing continuous deterministic functions (for the first 8, 
the fluid limits in the translation terms of (4.5)), with two null components 
§2,1 ee Oe and D = Oe. Five components of the limit are determined by the 

relations Qi = PiQ s , Z^ t i = = 0e and Z 2i 2 = —Zip. Finally, (Q s ,Zi >2 ) 
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is the unique solution of the following two-dimensional stochastic integral 
equation: 

Q s {t) = Q s (0) + (>2,2 - ^1,2) / £1,2(8) ds - (pi0i + p 2 2 ) / Qs(s)ds 

jo io 

(4.8) +Li(t) - Li )2 (i) - 5i, 2 (t) - L 2j2 (t) - S 2)2 (t), 

Z l}2 (t) = Zi, 2 (0)-/ [(^2,2 -^1,2)^1,2(^(5)) +^1,2] £1,2(5) ds 
Jo 

(4.9) -Li, 2 (i) + L2,2(t) + ^2W, 
where, for i = 1,2, 

L 1 = A 1 +A 2 -U 1 -U 2 - 5i,i = {B 1 (71 (t)) : t > 0}, 

(4.10) L- 2 = {J3 i)2 (&, 2 (*)) : f > 0}, Si, 2 = {B ii3 ( 7ij2 (t)) : i > 0}, 

£2 = {#2 (72W) : t > 0} and / = {B 2 ( 73 (0) : i > 0}, 

with B\, Bi >2 , B 2>2 , i?i,3 ; -62,3 and B 2 being six independent standard BM's, 
while 7j, 7i,2 and (f>i 2 are strictly increasing continuous deterministic func- 
tions. Specifically, 

71 (t) = (Ai + A 2 + m\nx,\)t + (pi^i + p 2 B 2 ) / q s {u) du 

Jo 

^1,2(0=^1,2 / (1 - K\ t2 {x(u)))zi t2 (u)du, 
Jo 

4>2,2{t) = /i 2)2 / vri, 2 (x(n))(m2 - «i,2(«)) du, 

(4.11) J ° 

71,2(0=^1,2 / Ki i2 (x{u))zx i2 (u)du 
Jo 

72,2(0=^2,2 / (1 -7ri, 2 (a;(n)))(m 2 - zi, 2 (u))du, 

JO 

72(0= / V ,2 (^0))'7 2 (2;('")) rf - ", 73(0= / 0- 2 {x{u))du, 

Jo Jo 

where 

(4.12) i>{x{u)) = n 2j2 (m 2 - zi >2 {u)) + (J,i l2 z li2 (u), u>0, 

with tti^ 2 {x{u)) and o~ 2 (x(u)) being the quantities associated with the FTSP 
D(x(u), ■), defined in (2.7) and (4.1), respectively, and characterized in 
(4.4). 
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Since the FCLT describes a refinement of the transient behavior of the 
fluid limit, it should not be surprising that the limiting stochastic process 
(Q s , ^1,2) would be difficult to analyze. On the positive side, we can solve 
for Z\ y 2 in (4.8) without having to simultaneously solve for Q s , but we need 
Z\ t 2 to solve for Q s . An additional complication for Q s is the dependence 
between the driving Brownian motions for the two processes Q s and Z\ 2\ 
note that the time-transformed Brownian terms appear in both. 

The FCLT shows the impact of system variability on the stochastic limit. 
First, and perhaps of greatest interest, there is a Brownian contribution 
L2 = B 2 (72(*)) from the FTSP appearing in the equation for Z\^\ n ote the 
dependence between L2 and /. However, (L2,/) is independent of all other 
Brownian terms. We thus see that the fluctuations about the fixed target 
ratio r in the queue-difference process (2.1) due to FQR do have an impact 
on the stochastic limit. 

On the other hand, we see that the stochastic fluctuations associated with 
external arrivals and abandonments only affect Q s ; they have no impact on 
Z\2- The same is true for the stochastic fluctuations of service facility 1, 
which is always busy, without any sharing. These fluctuations are captured 
by the Brownian term L\ = B\ (71 (t)). However, as noted above, in distinct 
contrast, the stochastic fluctuations in the service processes at service facilty 
2 have a more complicated impact, because they appear in the Brownian 
driving processes of both equations. 

4.3. Important Corollaries. The stochastic limit in the FCLT depends 
critically on the fluid limit x, which typically must be computed numeri- 
cally, but an efficient algorithm was developed in [23], exploiting the QBD 
structure of the FTSP D when n 2 is rational. Since we are mainly interested 
in the steady state variance of the diffusion limits, and since the stochas- 
tic fluctuations become more significant when the fluid is nearly constant 
(which happens when it is close to its stationary point) it is reasonable to 
initialize the fluid model at this fluid stationary point in order to simplify 
the expressions in (4.8) and (4.11). We do this in the next corollary. 

From an application point of view, the fluid limit is "more important" 
than the refined stochastic limit during the fluid transient period, since then 
the changes in the prelimit are of order 0(n). It follows from Theorem 3.2 
that after some (relatively short) time, the fluid stabilizes close to its unique 
stationary point Xq in (3.3). After that happens, the refined stochastic limits 
become the significant approximation to consider. 

When we consider the stochastic refinement of the stationary fluid limit 
x* , the stochastic limit process becomes much more tractable: it is a bivari- 
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ate Ornstein-Uhlenbeck (BOU) process centered at the origin, as in [2, 30]. 
Consequently, the random vector (Q s (t), Zi j2 (t)) has a bivariate normal dis- 
tribution with zero means for all t, and the associated steady-state random 
vector (Q s (oo), Zi t2 (oo)) can be very useful in applications. It is character- 
ized by three parameters: the two variances and the covariance, which we 
exhibit explicitly in (4.16) below. 

For a matrix M, let M t denote its transpose. The following is the key 
result for applications. It gives explicit Gaussian approximations for the 
steady-state distributions of all quantities of interest. 

Corollary 4.1. (FCLT with a stationary fluid) If, in addition to the 
conditions of Theorem £.1, x(0) = x* for the stationary point x* in (3.3) 
so that x is stationary, then the time transformations in (4.11) simplify by 
having n(t) = &i 3 % 2 (t) = £i,2*, and (j>i, 2 (t) = mrf, i = l,2, where 

£i = 2(Ai + A 2 ) - Mi,2«i,2 ~ V%2(m 2 - z* j2 ), 

6,2 = ^1,2^1,2 (^*)z*,2, 6,2 = ^2,2(1 ~ 7Ti j2 (x*)(m 2 ~ Z*^), 

(4.13) 171,2 = /^i,2(l - 7ri,2(a;*))^*,2> 
172,2 = ^2,2^1,2 (x*) (m 2 - z* >2 ), 

C2=ip 2 (x*)a 2 (x*) and & = a 2 (x*), 

for a 2 (x*) and ip(x*) defined in (4.1) and (4.12) with x(u) = x* . Then 
(Q s , £1,2) becomes a BOU process, satisfying the two-dimensional stochastic 
differential equation (sde) 

(4.14) dX = MX + SdB, 

where X = (Q s , Z\^f , B = (B\, B 2 ) , with B\ and B 2 being two independent 
standard BM's, and 

Mi,i = -(pi0\ +p 2 9 2 ), Mi, 2 = (n2,2 -1*1,2), M 2 ,i = 0, 

-^1,2^2,2^24,2 

M 2 ,2 = ; ; 7 —TT < 0. 

+ l*2,2(m 2 - Zl )2 ) 

(4.15) 5j x = 6 + 6,2 + 6,2 + 171,2 + m,2 = 2(Ai + A 2 ), 
«->i,2 = 5 2 ,i = 171,2 — 172,1 = 0, 5 22 = £ 2 + £4, 

2fJ,i, 2 n 2 , 2 zt 2 (m 2 - 4 2 ) 



6 = 171,2 + 172,2 



Mi,22i 2 + (m 2 - z\ 2)1*2,2 ' 



As a consequence, (Q s (t), Z\ )2 (t)) has a bivariate normal distribution with 
zero means for each t. The covariance matrix of the steady-state random 



18 



O. PERRY AND W. WHITT 



vector (Qs(oo), ^1,2(00)) has elements 
(4.16) 

a 2 Qs (oo) = Var(Q s ) = Q l + Q 2 , 

^1,1 / Ai + A 2 



Q 2 



^l^LzuW = / (M2,2-/il,2)^ a , Zl|2 (00) ' 
|A^l,l| I Pl#l+P2#2 



52 

z _ ^4 _ 4,2 _ g 2 _V^>V) 

1 2|M 2)2 | m 2 ' 2 2|.M2,2| " 2|7W 2 ,2| ' 

4.,z ll9 (°°) = Cot-(Q s ,Z li2 ) = ^^^(oo), £5 = ( \ Ml M + 2 M22 \ ) ■ 

Proof. By the definition of a stationary point, if x(0) = x* then x(t) = 
x* for all i > given in (3.3); then 7Ti 2OE*) appears in (3.4). The expressions 
in (4.13) follow directly from the expressions in (4.11), by replacing the time- 
dependent fluid quantities by their stationary counterparts. The resulting 
pair of integral equations for (Q s (t), Zi t2 (t)) is known to be equivalent to the 
BOU sde in (4.14). The covariance matrix of the stationary distribution, E, 
is known to satisfy the matrix equation ME + EM* = — V, where V = SS l , 
from which (4.16) follows; e.g., see [2, 16]. Algebra shows that £4/2|.A/f 2,2! = 
(1 - (zf t2 /m 2 )). ' □ 

Remark 4.1. (when components become null) Notice that the results in 
Corollary 4.1 simplify greatly with pool-dependent service rates, i.e., when 
Mi, 2 = 1^2,2 - Mi,2 = 0. Then Q 2 = and £5 = 0, so that 0q S) z 12 (°°) = °- 

We now see how Theorem 4.1 simplifies under the condition of pool- 
dependent service rates (no longer assuming that x(0) = x*). 

Corollary 4.2. (FCLT with pool-dependent service rates) If, in ad- 
dition to the assumptions of Theorem 4-1, ^2,2 = Mi,2 — v , then the two 
diffusion-limit processes Q s and Z\^ 2 can both be represented as separate 
one- dimensional processes, which satisfy the following integral equations 



(4.17) 



Qs{t) = 4(0) - m [ Q s (s)ds + B 1 (^ 1 (t)), 
Jo 

Zi, 2 (t) = Z li2 (0) - v I Z lj2 (s) ds + B 2 (7 2 (i)) , 
J 
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where 

71 (t) = 2(Ai + A 2 )f +{j 2 ~ 9«(0)J (1 " e -<M ) 
72(f) = U (J [ m 2^i,2{x{u)) +z\$(u) - 2ir 1)2 (x(u))z 1)2 (u)]dii\ + 72(f), 
rai/i 72(f) defined in (4.11) ; 

(4.18) 771 = Ai + A 2 - mi/Ji,! - m 2 z^, r/ 2 = + p 2 6> 2 , 

6ui £>i and B 2 are dependent standard BM's. 

Proof. It is immediate from the expression for Q s in (4.8) that when 
A*i,2 = ^2,2 the diffusion process Q s can be analyzed separately from Z± t2 . 
Since qi = piq s and /Ui j2 = f-2,2, it follows from (3.2) that q s (t) satisfies the 
simple ordinary differential equation 

q s (t) = (Ai + A 2 - mi/ii,! - m 2 // 2)2 ) - (pi#i + £>2#2)<7s(f) = »7i - mq s (t), 
whose solution is 



«.(*) = - + ?(o) - ^ 

for 7/1 and fj 2 in (4.18). Notice that 71(f) here corresponds to 71(f) +71, 2(f) + 
72,2 (f) in (4.11). Inserting q s (t) above into 71(f) in (4.11) gives 71(f). Notice 
that 72(f) corresponds to 0i,2(f)+</>2,2(f)+72(f) in (4.11). Again substituting 
yields the conclusion. □ 

Remark 4.2. (Equivalence with the single-class model.) If, in addition 
to the conditions of both Corollaries 4.1 and 4.2, we also have 0\ = 6 2 = 0, 
then the diffusion-limit process Q s is the same as the limit obtained for the 
M/M/n + M model in the efficiency-driven (ED) regime, see [33]. That is, 
Q s is an Ornstein-Uhlenbeck process with infinitesimal mean equal to 6 and 
infinitesimal variance 2A = 2(Ai + A 2 ). Thus, its steady-state distribution is 
normal with mean zero and variance X/9. However, Z\ )2 remains somewhat 
complicated involving 72(f) in (4.11). 

4.4. The Case r = 1: Longer Queue First (LQF). The most complicated 
feature in the FWLLN and FCLT asymptotic results in the previous two 
sections, inhibiting application, is the need to analyze the FTSP. Specifically, 
both the approximating fluid model and the stochastic refinement depend 
critically on the FTSP D = £(7) = {#(7, s) : s > 0} at each point 7 G A. 
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In particular, both limits depend on .0(7) through the two functions 7Ti 2(7) 
and cr 2 (7). These two functions can be computed numerically, as indicated 
above. For the stationary fluid point x*, tt\ 2^*) is given explicitly in (3.4). 

However, there is an important special case, itself of practical value, in 
which the analysis simplifies greatly, which can provide insight more gen- 
erally. When the target queue ratio is r = 1, the FTSP D{p/) becomes an 
ordinary birth- and- death (BD) process for each 7 6 A. Then the quantities 
711,2(7) and a 2 {^) are both easily expressed. It turns out that they can be 
expressed in terms of the first two moments of the busy-period distributions 
of two M/M/l queues. We consider that case now. 

We now assume that r = 1, and take 7 € A. In this case, the FTSP evolves 
as one BD process when .0(7) > and evolves as another BD process when 
D("f) < 0. We call the boundary state. Let Ai(7) denote the constant rate 
up (away from the boundary) and let ^1(7) denote the constant rate down 
(toward the boundary) of ^(7) when D(^) > 0. Focusing on the movement 
relative to the boundary, let A2(7) denote the constant rate down (away 
from the boundary) and let ^2(7) denote the constant rate up (toward the 
boundary) of ^(7) when D(j) < 0. 

Note that we need to analyze .0(7) only through the associated stochastic 
process 

(4.19) X( 7 ,t) = l {D(jit)>0} , t>0, 

which records which region D(j,t) is in at each time t. The stochastic pro- 
cess X = X(j) = {X(j,t) : t > 0} in (4.19) is a {0, l}-valued process 
associated with an alternating renewal process. Let ^1(7) denote a time in- 
terval between the instant of a state change from state to state 1 until 
the next instant of a state change from state 1 back to state 0. Similarly, let 
^2(7) denote a time interval between instant of a state change from state 1 
to state until the next instant of a state change from from state back to 
state 1. The successive times in the alternating renewal process are indepen- 
dent random variables distributed as ^(7) and ^2(7). The process X(j) is 
a regenerative process in which the regeneration times can be the successive 
instant of a state change from state to state 1 until the next instant of the 
same state change again at a later time. The intervals between successive 
regenerations are distributed as T\ (7) + T2 (7) . 

Now observe that ^(7) is distributed as a busy period in an M/M/l 
queue with arrival rate Aj(7) and service rate ^(7), i = 1,2. In this context, 
the condition 7 G A is equivalent to Aj(7) < ^(7), i = 1,2. Under this 
condition, ^(7) is known to have a finite moment generating function with 
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a positive radius of convergence, so that all moments of Tj( 7 ) are finite. Let 

(4.20) mi(7) = l/m(j) and Pi ( 7 ) = ^(7)^(7), i = 1,2. 
Then, from basic M/M/l theory, we have 

(4 . 21) mm-^ mm-^^. 

Finally, we are interested in the cumulative process associated with X(p/) 
in (4.19), 

C(l,t)=[ X(j,s)ds=[ 1{D( 7 , S )> } ds, t>0. 
Jo Jo 

We can apply (4.4) to obtain the following result. 

Theorem 4.2. (the FTSP when r = 1 ) When r = 1 and 7 6 A, the 
FTSP becomes a recurrent BD process. Hence the key FTSP quantities can 
be expressed directly in terms of the four BD rates Aj( 7 ) and ^(7) via 

(4-22) MT)^^«L, * 2 (7)- 



£[Ti( 7 )] + £[T 2 ( 7 )] w/ £?[Ti( 7 )] + S[T 2 ( 7 )] 
/or E[Ti(j)] and E[Ti(j) 2 ] in (4.21) and (4.20), i = 1,2. 

In the more general QBD setting arising with r 7^ 1, the analysis is more 
complicated, because the excursions of J * 1{_d( 7jS )>o} above and below de- 
pend on the entering and exit states from level 0; thus these excursions are 
not simply independent. Theorem 4.2 can be the basis for heuristic exten- 
sions to non-Markovian models in which the arrival, service and abandon- 
ment processes are non-Markovian. We may then exploit approximations for 
the busy period in GI/GI/1 queues, e.g., [1, 25]. 

5. Proof of Theorem 4.1. First observe that the assumed convergence 
in R2 a t time in (4.6) is actually equivalent to the full convergence in R17 
of the process in (4.7) at time because of Assumption 5. Our proof has 
four main steps: The first step is to exploit SSC results established in [24]. 
In particular, we first give an asymptotically equivalent three-dimensional 
representation of Xq (without any scaling) involving rate-1 Poisson pro- 
cesses. Then we observe that the essential dimension is actually two (when 
scaling by y/n) because the queue lengths are asymptotically in the fixed 
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ratio. Thus we deduce that it is sufficient to directly prove convergence of 
the 2-dimensional process {Q™,Z\2)- 

The second step is to facilitate application of the continuous mapping 
theorem by showing that an essential mapping is continuous. The third step 
is to construct appropriate martingale representations, allowing application 
of the continuous mapping theorem. The fourth and final hardest step is 
to show that the driving stochastic terms in this martingale representation 
converge to the specified limits. This final step uses a new result of inde- 
pendent interest, Theorem 6.1, the generalization of the classical FCLT for 
cumulative processes in (4.1) to the case where the QBD parameters at time 
t are given by the fluid limit x(t), which in general is time- varying. 

5.1. Representation and SSC. Following common practice, as reviewed 
in §2 of [20], we represent the processes A™(t), S™j(t) and U™(t) introduced 
at the beginning of §4.2 in terms of mutually independent rate-1 Poisson 
processes; let 

A?(t) = N?(\?t), 

(5.i) ^ = n * si z? ' j{s) ds ) and sn ^PSt S *' P 

U?(t) = N?(0ijQ?(8)d8Y t>0, 

where N*, Nf- and Nf for i = 1,2; j = 1,2 are eight mutually indepen- 
dent rate-1 Poisson processes. Theorem 5.1 of [24] gives a representation of 
the CTMC in terms of these processes. Corollaries 6.1-6.3 plus Theorem 6.4 
of [24] then establish state space collapse (SSC) results yielding an asymp- 
totically equivalent three-dimensional representation of Xq involving these 
mutually independent rate-1 Poisson processes plus two others. Since we ex- 
ploit that representation, we state it here. Directly, the representation of Z™ 2 
below keeps it in the interval [0, m^]. However, the representation directly 
allows the queue lengths Qf to become negative. The results in [24] show 
that the occurrence (anywhere in a a bounded interval) is asymptotically 
negligible. Recall that dj 1 denotes the Skorohod J\ metric. 

Lemma 5.1. {Representation via SSC of the service process) Under the 
assumptions in Theorem 3.1, dj x (Xg , Xq'*) =4> in Dq as n — > oo, with the 
three determining components of X^'* in (2.3), i.e., in X n in (2.4), being 
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represented via 



z i,2(t) = ^"2(0) + / 1 {D{\(s-)>o} dS 2j2 (s) ~ / l{D" 2 (s-)<o} dS^ 2 (s) 
Jo Jo 

± + Ni, 2 L 2 ,2 J 1 { D» 2 ( S )>0}(^ - ^2 W) ^ 



Q?(i) = Q?(0) + AJ(t)- / 1 W2(S _ )>0} ^( S )- / l W .-)<o}dS?,i(s)-E?(i) 







4 Q?(0) + JV?(A?t) - N^^t) - JVg U )2 jf l {D n 2(s)>0} ^2(s)) d« 

" #2,2 6*2.2 jf 1{D?, 2 ( S )>0}(™2 " ^I n 2(^)) df) " #1 (*1 jf Ql (*) &) , 

Q2 (*) = Q2 (0) + ^(t) - |* l{D^ 2 (,-)<0} dS^W - |* l { D ri2 ( s -)<o} - C/ 2 n (t) 

= QJ(0) + iV 2 a (A^) - iVg ^ 2 , 2 jf* l {o?2 ( s )< }(m^ - Z? )2 (a)) ds' 



iV 



1,2 



Ml,2 J* l{D l2 (s)<0}Z^ 2 ( S ) d S ) - iV 2 " (e 2 £ Q n 2 (s) d S ) . 



where N^ 2 and N 22 are two additional rate-1 Poisson processes, indepen- 
dent of the others. 

The representation in Lemma 5.1 provides important simplification, but it 
also shows the difficulty in proving heavy traffic limit theorems; the integrals 
contain the indicator functions depending on D™ 2 . We now show that the 
essential dimension can be reduced from three to two when we introduce 
scaling. The next result follows from Corollary 4.2 of [24]. 

Lemma 5.2. (SSC to two dimensions) Under the conditions of Theorem 
3.1, the essential dimension can be reduced from 3 established in Lemma 
5.1 to 2, because dj t (Q™, rQ 2 )/a n =4> in T>([0,5)) for 5 in Theorem 3.1 
whenever a n /\ogn — > 00 as n — > 00. If x G A over an interval \t\,t 2 ), 
< t\ < t 2 < 00, then the conclusion holds in D((ti,t 2 )). 

Due to Assumption 6 and Lemma 5.2, it is sufficient to directly prove 
convergence of the 2-dimensional process (Q™,Z± 2 ); the more general 16- 
dimensional limit in (4.7) can be obtained as a byproduct of the analysis, 

and in particular, Qi = piQ s , i = 1,2. 
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5.2. A Continuous Mapping. As in [20], our proof exploits the contin- 
uous mapping theorem. However, in our case, the stochastic processes de- 
scribing the evolution of the system (the queue length and service processes) 
cannot be expressed directly as a continuous mapping of the primitive pro- 
cesses. We next establish the continuity of the mapping that we will even- 
tually apply. 

Lemma 5.3. (Continuity of the two-dimensional integral representation) 
Consider the two-dimensional integral representation 

xi(t) = b\ + y\(t) + a 2 / x 2 (s) ds + a\ / x\(s)ds 
(5.2) /° h 

x 2 {t) = b 2 + y 2 (t) + / g(s)x 2 {s)ds 
Jo 

where g : R — > R satisfies c/(0) = and is Lipschitz continuous with a 
Lipschitz constant c g . The integral representation (5.2) has a unique solution 
(x\,x 2 ), so that the integral representation constitutes a function f : T> 2 x 
R 2 —?■ T> 2 mapping (yi,y 2 ,bi,b 2 ) into (x\,x 2 ) = f(yi, y 2 , h, b 2 )- In addition, 
the function f is a continuous mapping from D 2 x M 2 to T> 2 . Moreover, if y 2 
is continuous then x 2 is continuous. If both y\ and y 2 are continuous, then 
x\ is also continuous. 

Proof. By the conditions on the function g we have for all T > 

llffllr < 5(0) + \\g(u) - g(Q)\\ T < 5 (0) + c g T = c g T. 

Note that x 2 does not depend on x\, hence we can prove the lemma iter- 
atively by first showing that the function f 2 : T> x R mapping (y 2 ,b 2 ) into 
x 2 = f 2(112, b 2 ) is continuous, and then use this result to show that the func- 
tion /1 : T>2 x R mapping (yi,X2,b\) into x\ = fi{yi,X2,b\) is continuous. 

To show that / 2 is continuous we use Theorem 2.11 in [27] with h(x2(u),u) = 
g(u)x 2 (u). For that purpose, choose T > and let A be a homeomorphism 
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on [0, T] with strictly positive derivative A. Then, for every (pi,(p 2 € T> 

g{u)ipi{u) - g{X{u))(p 2 {X{u))\ du 
< / \g{u)(fi(u) - g(u)<p 2 (X(u))\ du + / \g{u)ip 2 {X{u)) - g(X(u))ip 2 (X(u))\ du 



<\\9\\t \(fi(u) - tp 2 (\(u))\ du + \\ip 2 \\T / \g(u) - g(X(u))\du 
Jo Jo 

< \\9\\t / \(fi(u) - f 2 (X(u))\du + CgT\\(p 2 \\ T \\X- e\\ T 
Jo 

= c\\\X - e\\r + c 2 I \(p\(u) — ip 2 {X{u))\du. 



where c\ = c g T\\(p 2 \\T and c 2 = \\g\\T- 

For x\ = fi(yi,x 2 ,bi) we can apply Theorem 4.1 in [20] with input 
y = yi + «2 Jq x 2 (u) du. It follows from Theorem 2.11 in [27] that if y2 
is continuous then so is x 2 . If, in addition, y\ is continuous, then y is con- 
tinuous and, by Theorem 4.1 in [20], so is x±. □ 

5.3. Martingale Representations. As in Theorem 6.3 of [24], we next ap- 
ply the representation in Lemmas 5.1 and 5.2 to obtain martingale repre- 
sentations for and Z™ 2 , but now we are interested in the FCLT instead 
of the FWLLN. We exploit martingale representations for the counting pro- 
cesses appearing in lemma 5.1 constructed from the rate-1 Poisson processes 
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Nf, N° 2 , N-' 2 2 and JVJ 4 , i = 1,2, in particular, 
(5.3) 

Mi^t) = N i,2 (^1,2 l{D^ 2 ( s )<o}- z 'r,2(' s )^ - 1^1,2 l {Di 2 (s)<o}Zi t2 (s)ds, 
M% 2 {t) = iV|, 2 L 2>2 jf l { ^ aW >o}(m5 - ^ 2 ( s )) cfe 

- ^2,2 f 1{D" 2 ( S )>0}(^2 - Z r,2(«)) ds, 
Ml,3(t) = N l£ (1*1,2 J 1 {Dl 2 (s)>0}Zi^( s ) d ^j ~ Ml,2 ^ 1{D" 2 ( S )>0}^™2('5)^ ; 

M£ 3 W = ^2,2 (m2,2 jf l { D r>2 ( s) <o}(m^ - Z£ 2 (*)) ds) 

- ^2,2 / l{D™ 2 (s)<0}( m 2 - ^™,2( s )) ds > 

J 

M^(i) = iV»(Ari)-A^, z = 1,2, 

M2.(t) = N?(0iJ Q7(s)ds S \ -GiJ Q?(8)da, i = 1,2. 

Lemma 5.4. (martingale representation for ) 
(5-4) 

Q?(t) = Q?(0) + o* 2)2 - Mi,2) / Zl 2 (s) ds - { Pl e l + P2 e 2 ) [ &(s) ds 

Jo Jo 

+ Ms(t) + op(l) as n->-oo, 

where M™ = M™/y/n for the martingale 

2222 
(5.5) M«(i) = £ M« (t) - £ M£(i) - £ M£(i) -J2 - M^(t). 

i=l i=l i=l i=l 

with respect to the natural filtration. 
PROOF. By Theorem 5.1, 

Qs(t) = Qs(0) + (A? + >%)t - m?A*i,i* - Mi,2 / ^ >2 («) ^ - /i 2 , 2 / Z% 2 (s) ds 

Jo Jo 

-01 [ QUs)ds-0 2 [ Q$(s)ds + M?(t), 
Jo Jo 
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for Mg(t) in (5.5). Observe that the indicator functions in the representation 
of X n in Lemma 5.1 do not appear in the representation of Q™(t). That 
simplifies the analysis. 

From (3.2) it follows that q s = qi + q2, the fluid counterpart of Q™, evolves 
according to the integral equation: 

q s (t) = q s (0) + (\x + \ 2 )t - m^mxt - zi^(u) du - /i 2 ,2 / z 2 , 2 (u)du 

Jo Jo 

— 01 qi(u)du — 02 \ q 2 (u)du, 



JO JO 

so that, substituting qi with piq s (u) and q 2 (u) with p2q s (u), 
q s (t) = q s (0) + (Ai + X 2 )t - /zi,imit - H 2>2 m 2 t 

+ (M2,2-Ml,2) / zi >2 (u) du - (pi0i +p 2 2 ) / q s (u)du. 
Jo Jo 

Then, by centering about nq s and dividing by y/n as in (4.5), we have 
qu^ = gn^ + [(A? + Ag ) - n(Ai + A 2 )]t Mi,i( m i ~ nm i)* 



n \/n 



Mi,2 Jo i z i, 2 ( s ) ~ nz iA s )) ds ^2,2 Jq {Z^ 2 (s) - nz 2 , 2 (s)) ds 



n \jn 
0i Jo (Qi( s ) ~ ™M S )) ds e * Jo (Q2 M " Mfgjf + 



By Assumption 1, the second and third terms in the expression above con- 
verge to zero. By Corollary 6.2 and Theorem 6.4 in [24], n~ 1//2 ||Z 2 2 — (m^ — 
Z™ 2 ) || =4> in 2? as n — > 00 so that z 2>2 = ?^2 — ^1,2- Also, (m 2 — nm 2 ) / v 7 ™ — ^ 
as n — )• 00 by Assumption 1. Hence, 



(5.6) 



Define 



Qs = Qs(0) + 0*2,2 - Ml,2) / ^2 W 

JO 

-^1 / QUs)ds-0 2 [ Q^{s)ds + M^{t) + o P {l). 
Jo Jo 



C s W = «(0) + fe-M f Zl 2 {s)ds-pi0i ! Q n s {s)ds 

Jo Jo 

-P202 [ &(s)ds + M2(t) 
Jo 
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By applying the SSC result in Lemma 5.2, we conclude that \\Q™ — QosIIt 
in T> as n — > oo for any T > 0. That completes the proof. □ 

We now turn to the process Z™ 2 . 

Lemma 5.5. (martingale representation for Z™ 2 ) 

ZUt) = ZlAO) ~ f [(M2,2 - Ml,2)7T 1|3 (x(fl)) + Ml,2] Zl 2 { S ) ds 



(5.7) 

+ L n + Mf + o(l), 
where L n = L n /^, Mg = M%/y/n, 

(5 g) i n (t) = f\hD U s)>v } ~ 7ri, 2 (x(s))]* n (s)) ds, 

* n (s) = M2,2(m^ - Z£ 2 (a)) + Ml,2^ 2 (^) 
and Mg is the martingale 

(5.9) M|(i) = Ml 2 {t) ~ M?4t). 

with respect to the natural filtration, where M 22 an d M™ 2 the martingales 
in (5.3). 



Proof. We start by rewriting the representation of Z™ 2 in Lemma 5.1 



as 



Zl,2(t) = ^I n 2 (0) " Ml,2 / (1 - 7T h2 (x( S )))Zl 2 ( S )ds 

(5.10) t Jo 

+ ti2,2 I 7r li2 (x(s))(m^-Z^ 2 (s))ds + L n + Ml 
Jo 

To achieve the diffusion-scaled process, we center (5.10) about nz\ j2 and 
divide by y/n, where, by (3.2), the fluid limit z\ % satisfies the equation 

zi,2{t) = 21,2(0) + A*2,2 / 7r 1)2 (a;(s))(m2 - z^ 2 (s))ds 
Jo 

-^1,2/ (1 -ni,2(x(s)))zi >2 (s)ds. 
Jo 

We get the representation (5.7) with the o(l) term replacing the deter- 
ministic term [(m 2 — nm 2 ) J Q * iri, 2 (x(s)) ds\/^fn < (m 2 — nm 2 )t/ 'y'ra, which 
converges to zero by Assumption 1. □ 
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5.4. Convergence of Stochastic Driving Terms. Given the representa- 
tions in Lemmas 5.4 and 5.5, we can complete the proof of the convergence of 
(Q™,Z™ 2 ) hi Theorem 4.1 by establishing convergence of the driving terms 
and applying the continuous mapping theorem with the mapping in Lemma 
5.3, i.e., with the following lemma, proved in the next section. We add an 
extra process, I n , also defined in (4.5), which is closely related to L n , but 
not directly needed to treat (<5™,^f 2 )- 

Lemma 5.6. (convergence of driving terms) Under the assumptions of 
Theorem 4.1, 

(5.11) (M?,M2,L n J n )^(M s ,M z ,L 2 J) in V 4 , 

where 

M s (t) = B x ( 7l (t)) - fl lj2 (7l, 2 (*)) " £2,2 (72,2(0) " Bl,3 (01, 2 (O) " #2,3 (02,2 
M Z (t) = #2,2 (02,2(0) " #1,2 (01,2 (*)), 

L 2 {t) = B 2 { l2 {t)) and I(t) = B 2 { l3 (t)), t > 0, 

B\, £?i,2; -82,2? -Bi,3j -B2.3 an d B 2 are independent standard BM's as in the 
statement of Theorem 4.1 and ji(t), ") 2 (t), 73(0, 71,2 72,2(0? 0i, 2 (O an <^ 
02,2(0 are ^ e increasing continuous functions defined in (4.11). 

5.5. Overall Proof of Theorem 4-1- We prove convergence of (Q n ,Z™ 2 ) 
by applying the continuous mapping theorem with the continuous function 
in Lemma 5.3, exploiting the representations in Lemmas 5.4 and 5.5 and 
the convergence established in Lemma 5.6. In applying Lemma 5.3, we rely 
heavily on Theorem 7.1 in [23], which establishes that tti <2 (-) is locally Lips- 
chitz continuous in A as a function of the fluid state x and is thus Lipschitz 
continuous over compact sets. Moreover, x(-) is itself Lipschitz continuous, 
as a function of the time argument s by Lemma 8.1 in [24]. It follows that 
tti, 2 (x(s)) is Lipschitz continuous as a function of the time argument s as 
well (using Assumption 6 implying that x lies entirely in A). Thus the proof 
of Theorem 4.1 is complete with the exception of the proof of Lemma 5.6. 
The next four sections are devoted to that proof. 

6. Proof of Lemma 5.6. This section is devoted to proving Lemma 
5.6. In §6.1 we apply standard arguments to establish the convergence of 
the first two martingale terms (M",M^). In preparation for treating the 
last two terms, in §6.2 we state two key results that we will use; they are 
proved in the following three sections. In §6.3 we establish convergence of 
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the last two terms (L n ,I n ). Finally, in §6.4 we establish joint convergence 
of all four terms by proving asymptotic independence of the last two terms 
from the first two terms. 

6.1. The First Two Terms in (5.11). We start by establishing conver- 
gence of the first two terms in Lemma 5.6, the two martingale terms. 

Lemma 6.1. There is joint convergence of the martingale processes 
(6.1) (M s n ,M|)=> (M S ,M Z ) in V 2 , 

where the processes are defined in (5.5), (5.9) and Lemma 5.6. 

Proof. Let 

M$(t) = (M^(t),M^ 2 (t),M£ 2 (t),M£3(t),M£ 3 (t)) in V 5 , 

Ml(t) = (M ai (t),M a2 (t)) and M£(t) = (M£(t),M™ (t)) in V 2 

for the martingale processes in (5.3). To compress the notation, for x G T>^ 
and t G [0,oo) fc , let x(t) = (x\(ti), x 2 {t2), ■ ■ ■ ,%k{tk))- We start by proving 
that 
(6-2) 

(M^,M%(t),M{}{t)) (BAiXfyBsMt^Bu (o J q(s)ds 
as n — > oo, where 

(6.3) 0i(t) = mifii,i(t), fa(t) = (/>i,2(t), <h{t) = Mt), 
Mt) = 71,2 (t), (f> B (t) = 72,2 (t), 

for and 7^2 defined in (4.11). Here BA(t), Bs(t) and Bjj{t) are vectors of 
independent standard Brownian motions. Using our compressed notation, we 
have At = (Ait, X 2 t) and 6q(s) = (0iqi(s),0 2 q 2 (s)). For example, S^(At) = 
(Ait), 5^2(^2*))) an d similarly for -Bs'(-) and Bjj(-). 
To prove (6.2), we apply the FCLT for Poisson processes. For the Poisson 
processes in Lemma 5.1, let 

K = N;{nt) - n \ M"^'"^"' and 



in Vq 



vn 
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Let M%(t), Mg(t) and Mjj(t) be the corresponding vector- valued processes. 
By the independence of all the unit-rate Poisson processes iVf(-), Nfj(-) and 
iV"(-), and the FCLT for a Poisson process, the following joint convergence 
holds: 
(6.5) 

(M^),Mg(i),M£(i)) (B A (t),B s (t),Bu(tj) in V 9 as n -> oo, 

where -Ba, an d By are, respectively, 2-dimensional, 5-dimensional and 
2-dimensional independent Brownian motions; see Theorem 4.2 and §9.1 in 
[20]. 

We now introduce random time changes. Let 
*A,i(t) = n 
*S )2 (t) = n 

*S, 3 (*) = « 
*S (4 (t) = n 

*s, 5 (*) = n 

(6.6) *^(t) = n 

By Assumption 1 on the arrival rates, <3?^ — > A«e in 2?, i = 1,2. From the 
initial conditions in the statement of Theorem 4.1, the fluid limit and the 
continuity of the integral mapping, it follows that =4> fa, 1 < i < 5, and 

^f/i(*) =^ ^» Jo 9*( s ) in £> as n — >■ oo. 

Let <fr^(t), and 3>jj(t) be the corresponding vector- valued processes. 

By Theorem 11.4.5 of [32], these limits hold jointly, yielding 

(6.7) (n(t), *§(*), *&(<)) =► (At, <f>(t), 9 jf g( s ) ds) in P 9 

as n — >• oo. By Theorem 11.4.5 of [32], the limits in (6.5) and (6.7) also hold 
jointly. By definition, 

(Mm,Mm,m{t)) = (i^(*x(<)),^(*s(t)),^(*&(*))) • 

Thus, the convergence in (6.2) follows from the continuity of the composition 
mapping at continuous limits, Theorem 13.2.1 in [32]. Finally, the conclu- 
sion of the lemma itself then follows from the definition of M™ and M% in 



"Vl,2 / l{D" 2 (s)<0} Z l,2( s )d s , 
JO 

"V2,2 / !{/?« ( s )>0}(™2 - -^r,2( s )) ds > 
JO 

'1*1,2 / l{D" 2 (s)>0}^r,2( s )^ S ) 

•» o 

_1 M2,2 I ^{Di 2 (s)<o}(m2 - Zl 2 (s))ds, 
Qf(s)ds, i,j = 1,2. 
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(5.5) and (5.9), and the continuity of addition under continuous limits, e.g., 
Corollary 12.7.1 in [32]. □ 

6.2. Key Supporting Results for the Last Two Terms. In §4.1 we indi- 
cated that the stochastic limit will depend on the FCLT for cumulative 
processes associated with the FTSP, as stated in (4.1). As indicated in §6 
of [23], the FTSP with fixed state 7 is a QBD; its parameters are given ex- 
plicitly in (13)-(17) of [24]. Since the FTSP D(j, •) is a QBD for each state 
7, it is a relatively simple regenerative stochastic process for each state 7, 
assuming that 7 makes the QBD positive recurrent. However, in our appli- 
cation, the fluid state is not fixed at 7, but is instead given by the fluid 
limit x(t), which is a function of time t. That means that the parameters of 
the FTSP are actually time-varying. By Assumption 6, the FTSP with fluid 
state x(t) is a positive recurrent QBD for all states x(t) considered. More- 
over, by Lemma A. 5 of [24], the infinitesimal generator and the asymptotic 
variance of the QBD are continuous functions of the underlying state x(t). 
Since the essential matrix structure (e.g., the dimension of the matrices) of 
the QBD's depends only on the rational ratio parameter 77^, and thus does 
not change, the QBD is characterized by only finitely many parameters. As 
a consequence, we can establish a variant of the FCLT in (4.1), allowing the 
FTSP to have a time-varying state. 

In our remaining proof of Lemma 5.6, in particular for Lemma 6.2 below, 
we will want to generalize the state of the QBD. The parameters of the QBD 
depend on the fluid state 7 = (q\, q2, z± t 2), but also on the rest of the QBD 
parameters, in particular, also upon ( = (\i,mj]i,j = 1,2). In order to 
establish Lemma 6.2 below, we will want to allow the parameters Aj and rrij 
to vary, because they vary with n in the many-server heavy-traffic scaling in 
Assumption 1. The QBD also depends on the other model parameters 9i and 
Hi j, but they are fixed, so we do not include them. Thus, we will consider 
the more general "full" parameter state function ij = (£, 7) for 7/ = n{t) and 
7 = j(t) above, which we understand to be an element of the functions space 
T>. We obtain a conventional stationary QBD model for each full parameter 
state rj(t). 

Now we will establish a FCLT for 

/•nt 

(6.8) C n {t; n) = n" 1 / 2 / (l {D ( 7?(s/n ), s ) >0 } - m^vis/n))) da, t > 0, 

J 

where the state function 77 is an element of T> and D(n(0),0) is some fixed 
finite initial value. The new FCLT is by no means immediate, because the 
function 77 makes the new process {D (77(5/71), s) : s > 0} appearing in the 
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integrand of (6.8) neither a QBD nor a regenerative process. Nevertheless, 
we establish the following generalization of the FCLT in (4.1). The proof is 
given in §8. 

Theorem 6.1. (FCLT for FTSP with time-varying parameter state) 
Consider the FTSP D as a function of its parameter state function n spec- 
ified above, where n is a function in T>. Suppose that the QBD D(rj(t), ■) is 
positive recurrent for all n(t), < t < T. Then 



with B being a standard BM and, for each u, a 2 (n(u)) is the asymptotic 
variance of the cumulative process with constant full parameter state ij(u), 
as in (4.1)-(4.4). 

For Lemma 6.2 below, we will also want to extend the FCLT in The- 
orem 6.1 to full parameter state functions that are suitably near a given 
deterministic one. For that purpose, we use the following elementary corol- 
lary to Theorem 6.1 and its proof. We use the Prohorov metric dp^O^i, Y2) 
characterizing convergence in distribution in T>([0, T]); see p. 77 of [32]. We 
say that a parameter-state function rj is positive recurrent if the associated 
FTSP D(rj, ■) is positive recurrent. 

Corollary 6.1. (continuity of the FCLT for the FTSP with time-varying 
parameter state) Consider the FTSP D as a function of its parameter state 
function n specified above, where the parameter state function n is a positive- 
recurrent element ofV. For all e > and T > 0, there exists 5 > such 
that, ifr]' is a parametric state function satisfying \\r] — t]'\\t < 6, then 7/ 
is positive recurrent for all t in [0, T] and d-p t T(C(-;rj),C(-,r]')) < e where 
C(-;n) is the limit process associated with D(n,-) defined in (6.10). 

In Lemma 6.2 below, we will apply Corollary 6.1 to random state func- 
tions fj n which converge weakly to n as n — > 00, i.e., for which fj n =^> 77 in T> as 
n — > 00. To do so, we need to connect the queue-difference processes D" 2 ap- 
pearing in I n in (4.5) to the FTSP. We do that via the associated frozen pro- 
cesses, introduced in §7.1 of [24]. The frozen process {D^ (X n (t) , s) : s > 0} 



(6.9) 



C n (-;n) C(-;rj) in V([0,T}) as n 00 



where C n (-;n) is given in (6.8) and 



(6.10) 
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corresponds to the queue-difference process D™ 2 starting at time t, condi- 
tioned on the state X n (t) at time t under the assumption that the transition 
rates are fixed or frozen at the rates associated with the initial state X n (t). 
A key property, for applying Theorem 6.1 and Corollary 6.1 above, is that 
the frozen process can be represented as the FTSP with modified parame- 
ters. To express the connection, we write the frozen process and the FTSP 
as functions of the parameters (Aj, rrij,"f, s). As in equation (48) of [24], we 
have the representation 
(6.11) 

{D](X?,m],X n (t),s) : s > 0} = {D(\f/n, m]/n, X n (t)/n, ns) : s > 0}, 

where Dj on the left of (6.11) is the frozen process described above, and D 
on the right of (6.11) is the FTSP. 

Like the queue-difference process, the frozen process has 0(n) transition 
rates, whereas the FTSP has O(l) transition rates, because of the time 
scaling in (2.5). Thus the time variable s on the right in (6.11) is scaled by 
n. 

However, we need to construct a process that is made up of different 
frozen processes over different subintervals. Thus, for each n > 1, we will 
construct a process that is a different frozen process over each successive 
interval of length 1/n, but identical to the queue-difference process at each 
interval endpoint. In particular, we will construct the overall frozen process 
by setting 

(6.12) D'f(t) = D n f {X n ({k - l)/n),i - (k - l)/n), ^— ^ < t < —, 

n m 

< t < T, where D™ is the frozen process defined above. That is, we use a 
different frozen state and thus frozen process for each interval [(k—1) /n, k/n) 
in [0, T]. As a consequence, the frozen process state for the process as a 
function of t is thus 

(6.13) Xf(t) = X n ([nt\/n), < t < T. 

As a consequence of (6.11)-(6.13) above, we can simply write 

(6.14) {D](t) :t>0} = {D(fj n (t), nt):t> 0}, 

where fj n is a random full parameter state function with the special param- 
eter function given in (6.11) above, with the frozen state at time t given 
by (6.13). Corollary 6.1 is relevant because, by virtue of Assumption 1 and 
Theorem 3.1, for each T > 0, we have 

(6.15) \\fj n — t]\\t => as n — > oo, 
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where r\ has fixed components Aj, mj and x(t), t > 0. 

Hence, the FCLT for fixed positive recurrent state function rj, which holds 
by Theorem 6.1, also holds with r\ replaced by fj n by virtue of Corollary 6.1. 
However, it remains to show that the newly constructed frozen processes 
approximate the queue-difference processes suitably well. For that, we will 
use a special coupling construction, similar to the coupling constructions 
used in [24]. The following result is proved in §7. 

Lemma 6.2. For each n, we can construct the new frozen processes 
defined by (6.11)-(6.14) on the same underlying probability space with the 
queue- difference processes D™ 2 so that, A" =4> in D as n — >■ oo, where 

(6.16) A n (t) = Vn~ [ (l{D» 2 ( s) >0} " 1 {m (s )>o}) ds > 1 > 0- 

J 

6.3. The Last Two Terms in (5.11). We now establish joint convergence 
of the last two terms in Lemma 5.6. 

Lemma 6.3. There is joint convergence of the last two terms in Lemma 
5.6, i.e., 

(6.17) (L n J n )^(L 2 J) in V 2 , 

where the converging processes L and I n are defined, respectively, in (5.8) 
and (4.5), while the vector limit process is (L2(t) , I (t)) = ( .E>2 (72 (i), -62(73^)) 
for L>2 a standard Brownian motion and (72 (i), 73 (i)) i> n (4-11), as in (4.10). 

Proof. We start by considering just I n . We make a change of variables 
in (4.5) to get 

1 f nt 

(6.18) I n (t) = ^J o [l{D? i2 (,/»)>0} " n^ixis/n))] ds, 0<t<T. 

From either the original representation of I n in (4.5) or the equivalent al- 
ternative expression in (6.18), the main line of the proof should be evident: 
We show that the time-scaled queue-difference process D™ 2 ( s / n ) i n (6.18) 
is asymptotically equivalent to the scaled FTSP D(x(s/n), s), making the 
expression in (6.18) be essentially of the form of C n in (6.8). If we could just 
directly make that substitution, then the desired limit I n I would be an 
immediate consequence of Theorem 6.1. However, the desired substitution 
is only valid asymptotically. We actually achieve the desired approximation 
by the FTSP indirectly by approximating the queue-difference process and 
applying Lemma 6.2 and Corollary 6.1 in addition to Theorem 6.1. 
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Having established the limit for I n , we turn to L n . From (5.8), we know 
that L n differs from I n by having the extra term \I' n in the integrand. How- 
ever, by the FWLLN, Theorem 3.1, \E ,n =4> tp as n — > oo, where ip(t) = 
1^2, 2(^2 — zi^it)) + ^1,2^1,2 (t). Hence, we first write 

(6.19) L?(t) = V^ [ [l{Di 2 (s)>0} ~ 7ri, 2 (x(a))]^(s) (is, t > 0. 

«» 

Since ip is continuous, we can approximate it uniformly closely by a piece- 
wise constant function, with all discontinuities occurring at multiples of a 
small positive e. Hence, by approximation, we can assume without loss of 
generality that 

[t/el+l 

(6.20) £?(*)= £ ^(t), 

3=1 

where ^ in (6.20) is a constant for each j and Ij'it) has the structure of I n 
over the subinterval [(j — l)e, je) and is outside that interval. Hence, we 
have convergence of L™, jointly with I n , by essentially the same argument 
as for I n . Finally, we can write 

(6.21) \\L n -Ll\\ T < ||* n -^||T||I n ||T^0 as n -»• 00, 

because ||* n - V'Ht => and I n => / as n — >■ 00, so that ||i n ||T => \\I\\t, 
where \\I\\t is finite by the continuous mapping theorem with the function 
|| • ||t- Hence, we obtain the claimed joint limit for (L n , I n ). □ 

6.4. Joint Convergence in Lemma 5.6. To complete the proof of Lemma 
5.6, it remains to show that the two limits established in Lemmas 6.1 and 
6.3 actually hold jointly. The two separate limits directly imply the associ- 
ated tightness, which in turn implies the tightness for the sequence of four- 
dimensional processes. Thus, to prove convergence it suffices to show that 
the limits of all converging subsequences coincide. We uniquely characterize 
the joint limit by showing that the limit for every convergent subsequence of 
the sequence {(M™, Mg, L n , I 71 )} must be of the form (M s , Mz, L, I), where 
(M s ,Mz) is independent of (L, I), having distributions as determined above. 
Thus it suffices to establish the following lemma. 

Lemma 6.4. (independent limits) The limits (M s ,Mz) and (L,I) for 
every convergent subsequence of the sequence {(M™ , M% , L , I n )} are inde- 
pendent. 
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In order to prove Lemma 6.4, we use the following lemma. 

Lemma 6.5. (basis for independent limits) Consider random vectors (I n , D n , V n ), 
where D n = f(V n ) for some function f, P(I n E B\V n = v) = P(I n £ 
B\D n = f(v)) for all Borel sets B almost surely with respect to dP(V n = v). 
If, in addition, D n =>• Oe, I n => I and V n =4> V as n — >■ oo, then I is inde- 
pendent of V . 

Proof. Let be a continuous bounded real- valued function on the ap- 
propriate space for i = 1,2,3. It suffices to show that E\g\(I)g 3 {V] = 
E[gi(I)]E[g 3 (V] for all continuous bounded real- valued functions g\ and 
$3. Since I n => I and D n => Oe, we also have (I n ,D n ) =^ (I, Oe) by Theorem 
11.4.5 of [32]. Thus, E[ gi (i n )g 2 (D n )} =► E[g x {i)g 2 {D)] = g 2 (0e)E[ gi (i)} as 
n — > oo. Then, by the assumptions above, using regular conditional proba- 
bilities since the stochastic processes are elements of the complete separable 
metric space T>, 

E[ gi (ng 2 (D n )g 3 (V n )] = J E[ gi (hg 2 (D n )g 3 (V n )\V n = v]dP(V n = v) 

= J E[gi(i n )g 2 (D n )\V n = v]g 3 (v)dP(V n = v) 

= J E[ gi (i n )g 2 (b n )\b n = f(v)]g 3 (v)dP(V n = v) 

= E[ gi (i n )g 2 (b n )]E[g 3 (V n )] 
-> E[ gi {i)g 2 {Qe)]E[g 3 {V)] 
(6.22) = g 2 (Oe)E\g 1 (i)]E[g 3 (V)] as n ^ oo. 

From the special case g 2 = 1, we see that E[gi(I n )g 3 (V n )] — > E[gi(I)]E[g 3 (V)] 
as n — > oo. Since the product gig 3 is a continuous bounded real-valued func- 
tion, we also have E[gi(I n )g 3 (V n )] — > E[gi(I)g 3 (V)] as n — > oo. Hence, 
E[gi(I)g 3 (V)] = E[gi(I)]E[g 3 (V)] for all continuous bounded real- valued 
functions g\ and g 3 , so that / is independent of V. □ 

Proof of Lemma 6.4. We show that the conditions of Lemma 6.5 are 
satisfied in our case. For that, we rely strongly on the SSC result in Theorem 
5.6 of [24]. We first observe that, for each n, the stochastic process {D™ 2 (t) : 
t > 0}, and thus also the stochastic processes {l{D™ 2 (i)>o} : ^ ^ 0} an d I n in 
(4.5) are directly functions of D n in (4.5). Thus, for each n, the conditional 
distribution of I n in V conditional on V n = (fo, UT 1 , S^,D n , Q™,Q n s , 
in (4.5) coincides with the conditional distribution of I n in T> conditional on 
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D n . Moreover, D n is scaled in the same way as the other processes in V n in 
(4.5). However, Theorem 4.4 (hi) of [24] implies that D n =>- Oe. Thus all the 
conditions of Lemma 6.5 are satisfied. Hence, that completes the proof. □ 

7. Proof of Lemma 6.2: A Coupling Construction. By the con- 
struction in (6.12), we have forced the new frozen processes to coincide 
with the queue-difference processes D" 2 for all time points t of the form k/n. 
To complete the proof, we employ a special coupling construction to con- 
struct these two processes on the same underlying probability space to make 
the processes have the same transitions within each interval [(k — l)/n, k/n) 
with high probability. As usual [19, 29], this coupling construction produces 
an artificial joint distribution, but leaves the distributions of each of the two 
processes individually unchanged. 

To complete the proof, we use the elementary bound ||A"||t < 2v n (T + 
e)/y/n for all n > 1/e, where T > and e > are arbitrary constants and 
where v n {T) is the number of k for which the nk <T and the sample paths 
of D r j and D™ 2 fail to be identical over the interval \{k — l)/n, k/n). We will 
show that the coupling construction can be done so that u n (T) = O(l) as 
n — > oo, thus implying the result. 

Since the transitions in the queue-difference process D™ 2 are generated 
from state changes in the CTMC Xq , we do the special construction from 
the perspective of the CTMC Xq. We use the device of uniformization to 
generate the transitions of the CTMC; i.e., we construct the transitions by 
thinning a Poisson process. Without loss of generality, we use different in- 
dependent Poisson process to generate potential transitions for each kind of 
transition, each interval [(k — l)/n,k/n) and each n. Since the transition 
rate of the CTMC is not uniformly bounded, there is a possibility that this 
direct construction will be invalid, but by choosing these Poisson process 
rates sufficiently high, we can make the likelihood of a violation asymp- 
totically negligible. In the actual construction, we can change the Poisson 
process when the constructed process hits a state from which a further tran- 
sition could lead to a violation. The detailed construction does not matter 
because we declare a difference occurring throughout the entire subinterval 
if the Poisson rate needs to be adjusted, thus contributing the maximum 
possible to the bound above. Since the integrand in (6.16) is bounded by 
2, the total impact upon (6.16) by such rate violations can clearly be made 
asymptotically negligible. 

The coupling is achieved by using the same Poisson processes to generate 
the transitions in both D" 2 and Df over each subinterval [(k — l)/n,k/n). 
These are done with respect to the states of X§(t) and Xg^k — l)/n). For 
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D™ 2 , the transitions rates of the various transitions (arrivals, abandonments 
from each queue and service completions of each class from each pool) are 
determined by the actual state X§(t), which changes throughout the interval 
[(k — l)/n,k/n). For, DJ, we do the same construction, but we leave the 
state fixed at its initial value X§((k — l)/n) throughout the interval \{k — 
1) /n, k/n), so that the transition rates do not change. However, we match the 
transitions in the two systems as much as possible. We make the transitions 
differ only to the extent that the state of X£(t) differs from Xg((k — l)/n). 

As stated above, we use different independent Poisson processes for each 
kind of transition. We have one Poisson process generate potential arrivals 
for each n. Since the arrival rates are unaffected by the state, the Poisson 
process for generating potential arrivals of class i can have rate Af , so that 
every potential arrival corresponds to an actual arrival in both systems. 
Thus no difference is caused by any arrival. That arrival in turn affects the 
constructed processes D™ 2 an d in the obvious way: an arrival of class 1 
increases them by 1, while an arrival of class 2 decreases them by r. 

For service completions of class 1 by pool 2, we let the Poisson process 
generating potential transitions have rate ui^m 2 . The actual transition rate 
at time t for X^t) is n\^Z™ 2 (t), so that the Poisson rate is an upper bound 
on the actual transition rate for all states. If the Poisson process with rate 
Mi, 2^2 h as a transition at time t, where (k — l)/n < t < k/n, then we 
let both systems have an actual service completion of class 1 by pool 2 at 
time t with probability {Z™ 2 {t) A Zf 2 ((A; — l)/n)/m 2 ; we let only the system 
associated with D™ 2 have an actual service completion of class 1 by pool 2 at 
time t with probability [Z^ 2 (t) - (Z^ 2 (t) A Z[\ 2 ((k - l)/ra)]/m£; we let only 
the system associated with DJ have an actual service completion of class 1 
by pool 2 at time t with probability [Z^ 2 ((k - l)/n) - (Zf 2 (t) A Zf 2 ((k - 
l)/n)]/m 2 ; and we let neither system have an actual service completion of 
class 1 by pool 2 with probability m£ - (Z^ 2 (t) V Zf 2 ((k - l)/n)]/m£. Thus, 
a difference in the sample path is caused by this transition with probability 
[{Zf 2 {t) V Z? 2 ((k - l)/n) - (Z^ 2 (t) A Z^ 2 ((k - l)/n)]/m£, which clearly is 
of order 0(l/n). 

We do similar constructions with independent Poisson processes for each 
of the other transitions. The abandonments are where the transition rate is 
unbounded, because the queue lengths Qf(t) are unbounded above. How- 
ever, for the Poisson process generating potential abandonments from the 
class-i queue over the interval [(k — l)/n, k/n), we can give it rate (Qf({k — 
l)/n) + 2Xi)0i, because the initial number in queue is Qf((k— l)/n) and new 
class-i arrivals occur at rate A". If the Poisson process generates a potential 
abandonment at time t, then it is a real abandonment for both systems with 
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probability [{Qf{(k - l)/n) A Q?(t))/(Q?((Jfe - l)/n) + 2A?); a difference 
in system states occurs only with probability [(Q"((fe — l)/n) V Q™(i)) — 
- l)/n) A Qi(t))]/(Qi((k - l)/n) + 2Af), which again is of order 

0(l/n). 

We now assess the impact of this construction. Both processes have tran- 
sition rates of order O(n) because the relevant processes and Zf- in Xq 
are O(n). Thus, the processes D™ 2 and have O(l) transitions over each 
interval of length 1/n. Hence, the state of X n (t) will only change an amount 
of order O(l) within each each interval [(k — l)/n,k/n). Consequently, the 
probability of any one transition being different is 0(l/n), and the proba- 
bility that there is any difference over the interval [(k — l)/n, k/n) is also of 
0(l/n). Hence, v n (T), the total number of intervals having any difference 
over the interval [0, T] will be of order 0(1), as claimed at the beginning of 
the proof. 

We can be more precise in this last step. We observe that, conditional 
upon Xq , which converges to xq, we can regard v n (T) as the sum of at 
most [nT\ + 1 independent Bernoulli random variable, assuming the value 
1 with probability p n ^ and otherwise, where Pi^/n < p n ^ < p u ,i/n for 
all i, provided that n is suitably large. The bounds hold because the set of 
values of (<?i, </2 5 -21,2) in x 6 fall into a compact subset of A and from the QBD 
continuity. Using the bound, we can bound v n (T) above stochastically by 
v™(T), defined as the partial sum of i.i.d. Bernoulli random variables, taking 
the value 1 with probability p u ,i/n. We can then apply the LLN for partial 
sums from triangular arrays to deduce that ^™(T) =>■ p u ,iT as n — > 00. That 
implies that is n (T) is indeed properly 0(1) as n — > 00. Hence the proof is 
complete. □ 

8. Proof of Theorem 6.1. First, if 7] is a constant function for which 
the QBD D(r],-) is positive recurrent, then the conclusion in Theorem 6.1 
reduces to the standard FCLT for a cumulative process in (4.1). To consider 
more general time-varying parameter state functions r\ = {i](u) : u £ /} 
(with I being the relevant time interval), we require that they be positive 
recurrent where, as before, we say that a state function r/ is positive recurrent 
if the associated QBD stochastic processes D(rj(u), •) are positive recurrent 
for all u £ I. Next, we observe that the conclusion in Theorem 6.1 is also 
valid for all positive-recurrent piecewise constant parameter state functions, 
where we include the condition that there be only finitely many disconti- 
nuities in each bounded interval. Let V pc be the subspace of V containing 
all such piecewise constant functions. To see that the conclusion holds for 
each positive recurrent rj G T> pc , note that, because of the time scaling, each 
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subinterval [a, b) of length 0(1) for the state function r/ corresponds to an 
interval of length 0(n) for the stochastic process {D(rj(s/n),s) : s > 0}, 
which has transition rates of 0(1). Moreover, the convergence on each suc- 
cessive interval implies that the initial distributions converge on the next 
interval. Hence, the initial conditions on each subinterval do not alter the 
limit. Thus, the separate subintervals can be treated separately, as if we 
were considering the first case of a constant parameter state function. 

Intuitively, it should be evident that the result extends to positive recur- 
rent state functions r] in D because each such function is the uniform limit 
over bounded intervals of piecewise-constant state functions; see p. 393 of 
[32]. However, a complete proof for this seemingly minor extension seems 
quite complicated. The remaining proof will be based on a series of lemmas, 
which are proved in the next section. 

First, we need to show that the subset of positive recurrent state functions 
in D is an open subset. In particular, we establish the following lemma. 

Lemma 8.1. (open subset) For any positive recurrent r\ in P([0,T]) ; 
there exists e such that if rf is another state function in D([0,T]) with 
II 7 ? — v'Wt < e » then rj is necessarily also positive recurrent in D([0,T]). 

With Lemma 8.1, we then exploit the continuity of QBD's established 
in Lemma A. 5 of [24] to complete the proof. We complete the proof in 
several steps, requiring further lemmas. First, as in the proof of Lemma 
6.2, we will exploit frozen processes to simplify the argument. As before, we 
use a coupling construction to show that they serve as suitable asymptotic 
approximations. Here we consider a modification of the process C n in (6.8), 
having a parameter state that is frozen over each successive cycle, where a 
cycle is the period between successive visits to a fixed state. In the case of a 
constant parameter state rj, these are the regeneration cycles associated with 
the regenerative process D(r/,-), as in [9], but here we have a more general 
case. For each n, let Cj denote this modification of C n . Closely paralleling 
the proof of Lemma 6.2, we use a coupling construction to show that it 
suffices to consider in order to establish the desired convergence of C n 
in (6.8). 

Lemma 8.2. (frozen cumulative processes) The processes C 1 ? and C n can 
be constructed on the same underlying space so that dj^C^j ,C n ) => 0. 

To prove CJ =4> C as n — > 00, we apply modified versions of the reasoning 
used to prove the FCLT in (4.1), as given in [9]. In particular, as in (1.1)- 
(1.4) of [9], we observe that C r f is asymptotically equivalent to a random 
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sum, ignoring remainder terms, and we then establish convergence for the 
sequence of random sums. To set the stage, let the i th full cycle in system 
n end at time T™. As usual, the key random variables associated with these 
cycles are the cycle lengths 

(8.1) T n = T n_ T n^ - > x 

and the integrals of the centered process over the cycle, which we call the 
cycle variables, 

(8-2) Yi 1 = / (l {D ( w , 8) >o} - 7ri, 2 (^)) ds, i > 0, 

where rfi = r)(T-]_ j), with TJLx being the random time at which the i th full 
cycle begins and Tff = 0, so that Yq is the cycle variable for the first partial 
cycle. We do not need to make additional assumptions for the analog of the 
variables Wi(f) in (1.2) of [9] because 

(8-3) W? = / |1 { d {??1 , s )>o} - 7ri, 2 (»7i) ds < 7?. 

With this construction, we can write 

(8.4) C]{T?;i 1 )=C n {T?;f l n f )i * > 0, 
for 

(8.5) *7?(<)=*fc, 27Li<t<3T. *<>o. 

Unlike in [9], here the random cycle vectors (t^Y™) are in general nei- 
ther independent nor identically distributed. However, the sequence of cycle 
variables {(r^YJ 1 ) : j > i} is conditionally independent of the entire sys- 
tem history up to time T^ 1: which we denote by J~lLi, given only T-^_ l7 for 
each i > and n > 1. Of course, in general these conditional distributions 
vary with i because the parameter state function r] is not constant, but they 
change little if 77 changes little, by the QBD continuity. Let N n (t) count the 
number of full cycles up to time t. As in (1.4) of [9], we can write 

(8.6) C]{t) = R n (t) + + K%(t), t>0, 
where R n (t) is the random sum 

N n (t) 

(8.7) R n (t) = n' 1 ' 2 ^ Y l\ t > 0, 

i=l 
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while Ri(t) and B%(t) are remainder terms involving the initial and final 
partial cycle, if any, also scaled by dividing by yfn. 

Just as in the standard regenerative setting, we are able to show that C" 
is asymptotically equivalent to R n , so that it suffices to work with R n . 

Lemma 8.3. [reduction to random sums) As n — > oo, i?" =>• Oe and 
R% Oe, so that dj^RJ^C]) => 0. 

It now suffices to show that R n (-;r]) =4> C{-\ rj) as n — > oo for each positive 
recurrent n in P. By virtue of Lemma 8.1, given such an rj, we can find a 
sequence of piecewise-constant state functions {r/™ : m > 1} where \\r]™ — 
ij\\t — 7-0 as m— t-oo with 77S being positive recurrent for all sufficiently large 
m. For those m, we have the desired convergence C n (-]rf p n c ) C{-;rj^) as 
n — > 00, as observed in the beginning of the proof. Thus, by Lemma 8.2 
and 8.3 above, we also have R n (-;ri™) => C{-',ff^) as n — > 00 for these m as 
well. We now want to show that the established convergence also holds when 
77™ is replaced by rj. For that purpose, we need to establish convergence as 
n — > 00 and m — > 00 jointly. In order to justify that joint convergence, we 
establish the following result. 

Lemma 8.4. (tightness and bounds for the random sums) Consider a 
parameter state function rj in T> and a piecewise-constant parameter state 
function rj pc , where both n and 7] pc are positive recurrent. Let T > and 
5 > be such that \\n — 7? pc ||r < S. Then the sequence {R n (-,r])} is C-tight 
in V([0,T*]) for some constant T* > and there exist functions ai(rj pc (-),5) 
and cr u (r)pc(-) , 5) such that the limit, say R(-,r]), of any convergent subse- 
quence of{R n (-,r])} can be represented as 

(8.8) R(t, V ) = B(W(t), V ), 0<t<T*, 

where B is standard BM and W can be bounded above and below by 

(8.9) f 2 a?( Vpc (s),5)ds < W(t2,ri)-W(ti,r]) < f \l{r] pc {s),8)ds 

for all t\ and ti with < t\ < ti < T* , where < af (rj pc (s) , 5) < 
o~u(r]p C (s),S) < 00 for all s, < s < T, and having the form in (4.4) deter- 
mined by the state rj pc (s). Moreover, for any e > and T* > 0, there exist 
5 > and T > as above, such that 

(8-10) \Wt(Vpc(-),6)-<ri(Vpc(-),6)\\T* < e. 
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Lemma 8.4 is based on associated lemmas for partial sums from trian- 
gular arrays of the cycle lengths and cycle variables r" and Y™, exploiting 
martingale structure; these results are stated in §9 and proved in §10. Given 
these lemmas, we now can complete the proof of Theorem 6.1. First, we 
have observed that C n (-,r] pc ) => C(-,r] pc ) in V for any positive-recurrent 
piecewise-constant parameter state function r/ pc . By Lemmas 8.2 and 8.3, 
R n (-,r]p C ) =^ C(-,r] pc ) in V as well. We can then apply Lemma 8.4 to de- 
duce that the sequence of random sums {R n (-,r])} is tight. Hence, each 
subsequence has a convergent subsequence. Let R(-,r]) be the limit of such 
a convergent subsequence. Next we construct a sequence {rj^} of positive- 
recurrent piecewise-constant state functions with t]\\t — > as m — > oo. 
As shown above, for each of them, we have {i? n (-,?/™)} =4> (7(-,ry™) as 
n — > oo. However, again by Lemma 8.4, we have R{-,rf) bounded above 
and below by the limits (7(-,?7™) which converge to {C(-,rj)} as m — > oo. 
Hence, we must have R(-,r]) = {C(-,r/)}. Hence all convergent subsequences 
must have the same limit, which implies that we must have the full conver- 
gence, R n (-,r]) C(-, rf) in T> as n — > oo. By Lemmas 8.2 and 8.3, we must 
also have C n {-,rf) =4> C(-,r/) in T>. Hence, Theorem 6.1 is proved. □ 

9. Proof of Lemma 8.4: Using the Martingale FCLT. We have 
indicated that Lemma 8.4 is based on associated lemmas for partial sums 
from triangular arrays of the cycle lengths and cycle variables r" and Y™, ex- 
ploiting martingale structure; in particular, we apply the martingale FCLT 
for triangular arrays. We can treat these two components of R (-,r]) sepa- 
rately because, just as in the familiar setting of renewal reward processes 
discussed in §§7.4 and 13.2 of [32], the FCLT for R n (-, rf) depends on a FCLT 
for partial sums of Y™ and a FWLLN for N n (t) separately. By the inverse 
relation discussed in §§7.3 and 13.6 of [32], a FWLLN for N n (t) is equivalent 
to a corresponding FWLLN for the partial sums of r". Since we can reduce 
the case of piecewise-constant rj pc to the case of constant r\ c by focusing on 
the subintervals separately, we now relate the given r\ to a constant r] c . 

Consider the cycle variables Yp in (8.2) associated with a parameter state 
function rj. Let J£ be the a-field generated by Xg(t) : < t < T™, k > -1. 
Let 

k 

(9.1) M^{k) = ^Y^,k> 1, and M$(t) = n~^ 2 M^([nt\ ), t > 0. 
i=i 
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For i > 0, let 

al^E[{Yrf\JF?-xl 

(9 2) L"*J 

V n (t) = n" 1 ^^ and V n (t) = sup{s : V n (s) <t}, t> 0. 
i=l 

We will be strongly exploiting the QBD continuity to obtain regularity in 
the variables Y™. 

Lemma 9.1. (sums of cycle variables) Consider a parameter state func- 
tion 7} in T> and an associated constant parameter state function rj c , where 
II 7 ? — Vc\\t < $ for some T > and 5 > 0, and both n and n c are positive 
recurrent. Consider the cycle variables Y™ in (8.2) and the associated vari- 
ables in (9.1) and (9.2), all associated with rj. Then there exist constants 
o~f(r] c , 5), o-^(r] c , 5) and 5' > such that, for all i and n, 

(9.3) o?(T) c ,S)<o* ti <ol(ri c ,S) and o*(t, c ,8) - af( Vc ,5) < 6' , 
for j in (9.2), associated with n, so that 

(9.4) afirjcS^-h) < V n (t 2 )-V n (t l )<al(r 1c ,5)(t 2 -t l ) 
for all n > 1 and < t\ < t 2 < T, for V n in (9.2). As a consequence, 

(9 .5) %zii) < v»(*a) - < {h ~ h) 



o-Kvab) tf(ilc,o~) 

for alln > 1 and < t\ < t 2 < T' for T' = T/a^(r], 5). Hence, the sequences 
{V n } and {V n } associated with n, defined in (9.2), are C -tight in V([0,T]) 
andD([0,T']), respectively. Moreover, the limits of convergent subsequences, 
say V and V must satisfy corresponding inequalities, i.e., 

af(Vc,S)(t 2 -h) < V(t2)-V(h)<al(ri c ,S)(t2-ti) and 

(9-6) i 4r z 4 Z V(fe)-V(t 1 )<%^ 

for the same ranges of t\ and t 2 above, so that V and V are both continuous 
and strictly increasing. In addition, 

(9.7) M$oV n ^B in V([0,T']) 



4G 



O. PERRY AND W. WHITT 



for My in (9-1), where B is standard Brownian motion. Thus, the sequence 
{My} is C-tight in T>([0,T]) with the limit of any convergent subsequence, 
say My, being of the form 

(9.8) M Y {t) = B(V(t)), < t < T, 

where V is bounded above and below over all subintervals as in (9.6). If we 
are free to choose the bounding constant 5 above, then for any e > 0, we can 
find 5 > so that 6' < e for 5' in (9.3). 

We now state the corresponding result for the partial sums of the cycle 
lengths. 

Lemma 9.2. (sums of cycle lengths) Consider a parameter state function 
rj in T> and a constant state function i] c , where both n and n c are positive 
recurrent. Let T > and 5 > be such that \\r] ~ Vc\\t < $■ Consider 
the cycle lengths r" in (8.1) associated with rj. Let U£ = rf + • • • + t^, 
k > 1, and U n (t) = n^U^,, t > 0. Let M™. = sfr^lJJLj], M"(t) = 

n~ 1 (M^ 1 H hM^ n(J ). Then the sequence {U n } is C-tight in V([0,T"]) 

for an appropriate time T" > 0, and if U is the limit of a convergent subse- 
quence, then necessarily it is bounded above and below with probability 1 by 
linear functions, i.e., 

(9.9) Pimiinc, 5)t < U{t) < m u {n c , 6)t, 0<t< T") = 1. 

where mi(rj c ,5) and m u (r] c ,5) are constants depending on 5 such that < 
< m u (n c ,5) < oo. If we are free to choose the time T > and the 
bounding constant 5 above, then for any e > and T" , < T" < oo, we can 
find 5 > so that the conclusions above hold with m u {r] c , 5) — mi{n c , 5) < e. 

As a consequence of the inverse relation between the partial sums and 
the associated counting processes, as in Chapter 13 of [32], we obtain the 
following corollary for the counting processes associated with the partial 
sums. Let N n (t) = n~ 1 N n (nt), t > 0. In the next section we combine 
Corollary 9.1 below with Lemma 9.1 to prove Lemma 8.4. 

Corollary 9.1. (counting process for cycle lengths) Under the assump- 
tions of Lemma 9.2, the sequence of scaled counting processes {iV n } is C- 
tight in D([0,T"']) for any time T'" < T"/mi(5), where T" is as in Lemma 
9.2. If N is the limit of a convergent subsequence of{N n }, then necessarily 
it is bounded above and below with probability 1 by linear functions, i.e., 

(9.10) P(t/m u (r) c , 6) < N{t) < t/m^ 6), 0<t< T'") = 1. 
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where mi(r] c , 5) and m u (j] c , 5) are the constants depending on 5 from Lemma 
9.2 above. If we are free to choose the time T > and the bounding constant 
5 above, then for any e > andT'" , we can find 5 > so that the conclusions 
above hold with m u (n c ,5) — mi(r) c ,5) < e. 

10. Remaining Proofs of Lemmas. In this section we prove the six 
lemmas in the previous two sections, which were used in the proof of Theo- 
rem 6.1. We prove them in the order needed for the proof. 

Proof OF Lemma 8.1. We exploit the criterion for recurrence in terms 
of the drift rates given in (2.8). The drift rates 5+(r]) and S-(jf) for constant r\ 
in the regions {s : D(rj, s) > 0} and {s : D(rj, s) < 0}, respectively, are linear 
functions of the components of the vector r\. We can thus express the drifts 
as the inner products S±(n) = a±- 77, where a+ and a_ are constant vectors. 
Hence, if \rj— rf\ < e, then | 6±(t])— S±(rf) | < e(|a±|T), where here 1 is a vector 
of l's of the appropriate dimension. This property for constant parameter 
states extends immediately to more general state functions in D using the 
uniform norm; i.e., if \\rj — t/\\t < e, then ||5±(?/) — S±(t]')\\t < e(|o±| • 1). 
Thus, for any positive recurrent state function rj, there exists e > such 
that 5+(r]') < and S-(rj') > if 1 1 77 — t]\\t < e, implying that 7/ is also 
positive recurrent. □ 

Proof of Lemma 9.1. The key observation is that the sequence of ran- 
dom vectors {(r™, : j > i} associated with the general parametric state 
function rj is conditionally independent of the entire system history up to 
time T™_^ for each i, which we have denoted by given only T-l_ 1 . As a 

consequence, paralleling the regenerative case in [9] and (4.4), 



(10.1) E 



' 1 {D(r,,, S )>0}) ^17^1 



for i > 1, where rji = r/(T i "_ 1 ), so that E\Y7 l \JFf_-^\ = for each i. Hence, the 
stochastic process {My(fe) : k > 1} is a square integrable martingale with 
respect to the filtration {J 7 ]} : k > 1}. 

Moreover, by the QBD continuity, the variances = E[{Y[ L ) 2 \ J-f_i\ 
in (9.2) cannot differ too much from the corresponding variance for the 
constant parameter function n c . For a fixed t > 0, let <Ty(7/(t)) be o\ i under 
the condition that T™_ x = t, so that rji = T)(T£_i) = r/(t). Since \\r]— T] c \\t < 5, 
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we can apply the QBD continuity to obtain the relations in (9.3), where 

vfiVciS) = min {oy (??(£)) : r/ £ A(r/ C ,5)} and 
a l(Vc$) = max {a Y (r)(t)) : rj G A(r) c , 5)} with 

(10.2) Afa c ,$) = {v:\\v-Vc\\t<6}. 

These in turn imply that the inequalities in (9.4) and (9.5) hold for V n and 
V™ for all n, implying the tightness of the sequences {V n } and {V n } and 
the inequalities stated in (9.6) for the limits of all convergent subsequences. 
However, we cannot conclude that in general either V n or V n converges. 

Nevertheless, we can apply an appropriate martingale FCLT to deduce 
that the limit in (9.7) holds; e.g., sec Theorems 2.1 and 2.2 of [5], Theorem 5 
of [26] and p. 98 of [14]. The QBD continuity and the bounds in (9.3) imply 
that the technical regularity conditions are satisfied in this case. Hence, for 
any 5 > 0, we can apply the martingale FCLT to get the convergence in 
(9.7). 

Given that V is a strictly increasing continuous function with bounded 
slope, as in (9.6), we can deduce from the tightness of {My ° V n }, which 
follows from the convergence in (9.7), that the sequence {My} itself must 
be tight. That is most easily done by letting V™ be a continuous function 
constructed by linear interpolation under which we still have the convergence 
in (9.7). Then, V n itself is a continuous strictly increasing function with 
modulus bounds in (9.6). Hence, we can deduce that the sequence {My} 
must be tight. 

The sequence {(My ,V n ,V n )} is tight because the component sequences 
are all tight. Starting from the joint convergence (My, V n , V n ) => (My , V, V") 
in V3 for any convergent subsequence, we can deduce from (9.7) that My = 
BoV, as claimed in (9.8). The final e bound follows from the QBD continuity 
in Lemma A. 5 of [24]. □ 

Proof of Lemma 9.2. The proof is similar to the proof of Lemma 9.1 
above, but now we need a FWLLN instead of a FCLT. However, it is con- 
venient to apply the FCLT in order to deduce the FWLLN. Indeed, by the 
same reasoning used to prove Lemma 9.1 above, we can obtain a martingale 
FCLT for the sums of the centered variables r™ — £'[r"|J r ^ 1 ], paralleling 
(9.7). Here we use the conditional variances and their sums, defined by 

[nt\ 

(10.3) 4, = EM ~ £[rf |JtiD 2 |^iL V n {t) = n- 1 J>^. 

i=l 



instead of (9.2). We then obtain analogs of (9.3), (9.4) and (9.6). 
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Given that FCLT, we scale further, essentially dividing by y/n, to get the 
associated FWLLN for the centered variables. As a consequence, we obtain 
the FWLLN U n - M% => Oe in V([0,T"]) as n ->■ oo, for an appropriate 
finite time T", not necessarily equal to T or T" in the previous proof above. 
Then, in direct analogy with (10.2), we apply the QBD continuity to obtain 
mi(Vc,8) < Mjjj < m u (rj c ,S) for all i and n. Hence, mi(r] c ,6)t < Mjj(t) < 
nT-uiVctyt f° r a h n an d t, < t < T. We can then combine these bounds 
with the FWLLN to obtain the conclusions stated in the lemma. By the 
QBD continuity, m u (r] c , 5) — mi(r] c , 5) — > as 5 ! 0. □ 

Proof of Lemma 8.4. First, Lemmas 9.1 and 9.2 and Corollary 9.1 can 
be extended directly to piecewise-constant state functions as well as constant 
state functions. Thus, for rj in V, they imply that the sequences {My} 
and {N n } are each C-tight in V. Consequently, the associated sequence of 
vector processes {(My , A? n )} is C-tight in T>2- Hence, every subsequence has 
a further convergent subsequence. Moreover, by Lemma 9.1 and Corollary 
9.1, any limit, say (My,iV), can be represented as {B o V,N), where V 
and N are bounded as in (9.6) and (9.10) over each subinterval where the 
piecewise-constant parametric state function is constant. Hence, overall they 
can be bounded above and below by 

(10.4) {V Y , h N t ) < (V Y , N) < (V Y ,u, N u ), 

where (^(0)^(0)) = (Vy, u (0), N u (0)) = (0,0) and 

(V Y ^(t),Ni(t)) = (V Y j(U-i) + a \i {P ~ ti-i), Ni(ti-i) + (l/m Uii )(i - U-i)), 
(V Y , u (t),N u (t)) = {V Y ,u(U-i) + <i(t - U-i),N u (ti-i) + (l/m M )(i - U-i)), 

for ti-i < t < ti, where = to < ti < • • • < tk = T, so that U are the 
endpoints of a piecewise constant state function rj pc , with af { and l/m u> i 
being the lower bounds and a^ li and being the upper bounds on 

the i th subinterval, depending on r] pc and 5. Hence, we can apply the con- 
tinuous mapping theorem to obtain the corresponding convergence for the 
random sum for all convergent subsequences, with the limit of all convergent 
subsequences represented as claimed in (8.8) with W there bounded as in 
(8.9). The bounding variance functions are given explicitly by &u(Vpc(s), 5) = 
a u,i/ m i,i an d (Vpc(s) , 5) = af Jm Ut i for U-i < s < U. Thus, by having 
II 7 ? — ^pcIIt < 8 and choosing S sufficiently small, we can obtain the desired 
variance inequality (8.10). □ 

Proof of Lemma 8.3. The reasoning follows the regenerative case as in 
[9]. First, the remainder term Ri(t) is relatively easy to treat since it involves 
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the initial cycle and is thus independent of t. Since D(rj(0),0) has been 
specified as some fixed state after (6.8), the initial partial cycle until hitting 
time of the designated state is clearly O(l) and becomes asymptotically 
negligible when we divide by \/n. 

As in [9], to treat the second remainder term, we exploit the representation 

< n- l / 2 W^ {t)+l <n~ l l 2 T n Nn{t)+1 
(10.5) < n _1 / 2 max{rf : 1 < i < N n (t) + 1}, t > 0, 

for W™ and r" defined in (8.3) and (8.1). However, the last term in (10.5) is 
asymptotically negligible because of the FCLT for the cycle lengths used in 
the proof of Lemma 9.2 above. The last term is the maximum discontinuity 
in the prelimit process indexed by n. Since the limit is continuous, that term 
is asymptotically negligible. □ 

Proof of Lemma 8.2. The reasoning here is similar to the proof of 
Lemma 6.2. As before, we can use a coupling construction to make the 
two processes have identical sample paths over the vast majority of the cy- 
cles. We exploit the oscillation property for functions in T>([0,T]), Corollary 
12.2 of [32], concluding that, for any e > 0, there are k time points t\ with 
= to < t\ < ■ ■ ■ < tfc-i < tk = T such that \r/(si) — 77(32) | < e for all 
si,S2 £ [ti_i,ti) for all i. Hence, with the time scaling by 1/n in (6.8), we 
see that, except for at most k cycles in [0, T] containing the k boundary 
points U, the oscillation of r\ over the cycle is at most e/n. Hence, the cou- 
pling can be performed as in the proof of Lemma 6.2, making the probability 
that the sample paths differ over any one cycle among all except the k be 
of order 0(l/n). Since there are 0(n) cycles in [0,T], as substantiated by 
Corollary 9.1, there are order 0(1) among the 0{n) cycles that have any 
difference in the sample paths. Hence, with the spatial scaling by y/n, we 
clearly have dj 1 (CJ,C n ) =4> as n — > 00 as claimed. □ 

11. Comparisons with Simulation. To both support the validity of 
the theorems and their applicability to the intended engineering problems, 
we now compare the approximations stemming from the FWLLN and the 
FCLT to the results of simulation experiments. Specifically, we will com- 
pare the Gaussian approximations for the steady-state queue lengths with 
simulation estimates of these quantities, obtained by simulating the actual 
queueing model over a large time interval. The approximate mean values 
come directly from the stationary point of the fluid limit, x* in Theorem 
3.2; the approximate variances come from Corollary 4.1, specifically, from 
(4.16). 
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Our simulation examples will have parameters related to a base case. 
First, scale is described by the parameter n, which is the scaling parameter 
in our limit theorems. The abandonment and service rate parameters, which 
describe the behavior of individual customers and servers, are independent 
of n: 9\ = 9 2 = 0.2, //^i = ^2,2 = 1-0 and p\ :2 = ^2,1 — 0.8. The service 
rates are chosen so that it is less efficient to serve a customer from a different 
class. 

The parameters that scale as the service system grows depend on n; they 
are chosen to be directly proportional to n: = nrrii, A- n ^ = n\i and 

k^l = wfei 2- We take k™ 2 to be order 0(n) so it is easy to compare different 
system sizes. Our base case then has mi = m 2 = 1, Ai = 1.3, A2 = 0.9 and 
k\ t 2 = 0.1. The arrival rates are chosen to put class 1 in a focused overload, 
while class 2 is initially normally loaded or slightly underloaded, but becomes 
overloaded too after the sharing. (These model parameters satisfy case 1 of 
Assumption 1 of [24].) We use the FQR-T control with ratio parameter 
r = 1.0, which allows us to apply the simple asymptotic formulas from §4.4. 

From (3.3) and (3.4), we see that the stationary fluid solution for this base 
case yields z\ 2 = 0.2111, q* = 0.6556, q 2 = 0.5556 and n* 2 = ni i2 (x*) = 
0.1763. Without any sharing, the fluid approximation for queue 1 would be 
1.5000. Hence the sharing reduces the first fluid queue from 1.5000 to 0.6556, 
at the expense of causing the second class to have a fluid queue of 0.5556. 

We now turn to the variances, for which we need to analyze the FTSP 
more carefully. The FTSP has BD parameters: Xi(x*) = 1.411, p,\{x*) = 
2.989, X 2 (x*) = 2.031 and fJL 2 {x*) = 2.369. The associated M/M/l traffic 
intensities are pi(x*) = 0.472 and P2(%*) = 0.8574. The associated mean 
busy periods are E[Ti(x*)} = 0.6338 and E[T 2 {x*)] = 2.9603. Hence, the 
alternative formula for ^1^(3;*) in (4.22) agrees with the value 0.1763 given 
above (providing a check on our calculations). 

Turning to the FCLT, from (4.12), we see that tp(x*) = 0.6200, so that 
^(x*) = 0.3844. For a 2 (x*), from (4.21), we see that E[Ti{x*) 2 } = 1.5218, 
so that Var{Ti{x*)) = 1.1201, and a 2 {x*) = 1.1201/3.5941 = 0.3116. Then 
6 = ^ 2 (x*)a 2 (x*) = 0.1198. Since \M 2)2 \ = 0.176, Z 2 = 0.3403. Hence, 
<r| ia (oo) = 1 - 0.2111 + 0.3403 = 1.1292.' 

As a consequence, <Tq SjZi2 (°°) = (1-1292)(0.5319) = 0.6006. Since p 2 ,2 — 
pi,2 = Pi0i+p 2 2 = 0.2, Q 2 = <Tq s ,z 12 (°°) = 0.6006. Since Qi = 11.0, we 
have <7q s (oo) = 11.6006, so that the associated standard deviation is 3.41. 
(Without Q 2 , we would approximate the standard deviation by \/TT = 3.32, 
so Q 2 contributes only 3% to the standard deviation approximation in this 
case.) 
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By the SSC, the diffusion approximations for Q\ and Q2 are linearly 
related to Q s ; in particular, Oq, (00) = {pi) 2 a g a (oo), so that Oq. (00) = 
11.6006/4 = 2.900 and the associated standard deviation is 1.70. 

We now turn to the simulations. We simulate the actual queueing system 
obtained by scaling up the appropriate parameters by n. We consider three 
cases: n = 25, n = 100, and n = 400. (Since k^ 2 must be an integer, we let 
fc" 2 = 3 when n = 25.) 

In all our simulation experiments, we used 5 independent runs, each with 
300, 000 arrivals. We report averages together with the half widths of the 
95% confidence intervals, based on a t statistic with four degrees of freedom. 
Simulation results for the base case above are presented in Table 1 below. 

The first four rows of Table 1 show mean values. We display both the 
steady-state mean values and the associated scaled values (i.e., divided by 
n). The unsealed values helps us evaluate the performance of the actual 
system, while the scaled values show the convergence in the FWLLN. Table 
1 clearly shows that the accuracy improves as n gets larger, but even for 
relatively small systems, the fluid approximation gives reasonable results. 





n=25 


n=100 


n=400 


perf. meas. 


Approx. 


Sim. 


Approx. 


Sim. 


Approx. 


Sim. 


E[Qi] 


16.6 


15.7 
±0.3 


65.6 


63.6 
±1.9 


262.2 


258.3 
±5.0 


E[Qi/n] 


0.656 


0.629 
±0.013 


0.656 


0.636 
±0.019 


0.656 


0.646 
±0.013 


E[Q 2 ] 


13.6 


15.9 
±0.4 


55.6 


58.6 
±1.8 


222.2 


223.9 
±5.0 


E[Q 2 /n] 


0.556 


0.636 
±0.016 


0.556 


0.586 
±0.018 


0.556 


0.560 
±0.013 


std(Qs) 
std(Q s ) 


17.1 


16.0 
±0.3 


34.1 


33.7 
±1.4 


68.2 


67.6 
±2.9 


3.41 


3.21 


3.41 


3.37 


3.41 


3.38 


std(Qi) 
std(Qi) 


8.5 


8.8 
±0.1 


17.0 


17.2 
±0.7 


34.0 


33.9 
±1.4 


1.70 


1.75 


1.70 


1.72 


1.70 


1.70 


std(Q 2 ) 
std(Q 2 ) 


8.5 


8.6 
±0.1 


17.0 


17.1 
±0.7 


34.0 


33.9 
±1.5 


1.70 


1.73 


1.70 


1.71 


1.70 


1.69 



Table 1 

A comparison of approximations to simulation results for the means and standard 
deviations of the steady-state queue lengths as a function of the scale parameter n in the 
base case with A" = 1.3n, \ 2 = 0.9n, fc" 2 = O.ln, r = 1 and other parameters defined 

above. 
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Rows 5 — 10 of Table 1 show the standard-deviations of the total queue 
length Q s = Qi + Q2 as well as the two queues. As before, we treat both 
the actual values and the scaled values, but now we are scaling in diffusion 
scale (dividing by \fn after subtracting the order-O(n) mean), as in (4.5), 
so that we will be substantiating the FCLT, specifically Corollary 4.1 and 
the variance formulas in (4.16). To save space, we omit the confidence in- 
tervals for the scaled standard deviations; these can be computed from the 
confidence intervals of the actual queues by dividing the half widths by \fn. 

Overall, we conclude that Table 1 shows that the approximations are 
remarkably accurate. 

Acknowledgments. This research began while the first author was 
completing his Ph.D. in the Department of Industrial Engineering and Op- 
erations Research at Columbia University and was completed while he held 
a postdoctoral fellowship at C.W.I, in Amsterdam. This research was partly 
supported by NSF grants DMI-0457095 and CMMI 0948190. 

References. 

[I] Abate, J. and Whitt, W. 1995. Limits and approximations for the busy-period 
distribution in single-server queues. Prob. Engr. Inf. Sci. 9 581-602. 

[2] Arnold, L. 1974. Stochastic Differential Equations: Theory and Applications, Wiley, 
New York. 

[3] Asmussen, S. 2003. Applied probability and Queues, second ed., Wiley, New York. 
[4] Coffman, E. G., Puhalskii, A. A. and Reiman, M. I. 1995. Polling systems with 

zero switchover times: a heavy-traffic averaging principle. Annals of Applied Probability 

5 681-719. 

[5] Durrett, R. and Resnick, S. I. 1978. Functional limit theorems for dependent 

random variables. Annals of Probability 6 (5) 829-846. 
[6] Fricker, C., Robert, P. and Tibi, D. 2003. A degenerate central limit theorem for 

single resource loss systems Annals of Applied Probability 13 (2) 561-575. 
[7] Gans, N., Koole, G. and Mandelbaum, A. 2003. Telephone call centers: tutorial, 

review and research prospects. Manuf. Serv. Oper. Mgmnt. 5 79-141. 
[8] Garnet, O., Mandelbaum, A. and Reiman, M. 2002. Designing a call center with 

impatient customers. Manuf. Serv. Oper. Mgmnt. 4 (3) 208-227. 
[9] Glynn, P. W. and Whitt, W. 1993. Limit theorems for cumulative processes. 

Stochastic Processes and Their Applications 47 299-314. 
[10] Gurvich, I. and Whitt, W. 2009a. Scheduling flexible servers with convex delay 

costs in many-server service systems. Manuf. Serv. Oper. Mgmnt. 11 237-253. 

[II] Gurvich, I. and Whitt, W. 2009b. Queue-and-idleness-ratio controls in many-server 
service systems. Math. Oper. Res. 34 363-396. 

[12] Gurvich, I. and Whitt, W. 2010. Service-level differentiation in many-server service 
systems via queue-ratio routing. Oper. Res. 58 316-328. 

[13] Halfin, S. and Whitt, W. 1981. Heavy-traffic limits for queues with many expo- 
nential servers. Oper. Res. 29 (3) 567-588. 

[14] Hall, P. and Heyde, C. C. 1980. Martingale Limit Theory and its Applications. 
Academic Press, New York. 



54 



O. PERRY AND W. WHITT 



[15] Hunt, P.J. and Kurtz, T. G. 1994. Large loss networks. Stochastic Processes and 

their Applications 53 363-378. 
[16] Karlin, S., and Taylor H. M. 1981. A Second Course in Stochastic Processes. 

Academic Press, New York. 
[17] Karr, A.F. 1975. Weak Convergence of a Sequence of Markov Chains. Z. 

Wahrscheinlichkeitstheorie und Verw. Gebiete 33 41-48. 
[18] LATOUCHE, G. and Ramaswami, V. 1999. Introduction to Matrix Analytic Methods 

in Stochastic Modeling. SIAM and ASA, Philadelphia. 
[19] Lindvall, T. (1992). Lectures on the Coupling Method, Wiley, New York. 
[20] Pang, G., Talreja, R. and Whitt, W. 2007. Martingale proofs of many-server 

heavy-traffic limits for Markovian queues. Probability Surveys. 4, 193-267. 
[21] Perry, O. and Whitt, W. 2009a. Responding to unexpected overloads in large-scale 

service systems. Management Sci., 55 (8), 1353-1367. 
[22] Perry, O. and Whitt, W. 2009b. A fluid approximation for service systems 

responding to unexpected overloads. Operations Res., forthcoming. Available at: 

http:/ /www. columbia.edu/~ww2040/allpapers. html 
[23] Perry, O. and Whitt, W. 2011. An ODE for an overloaded X model involving a 

stochastic averaging principle. Stochastic Systems, vol. 1, 17-66. 
[24] Perry, O. and Whitt, W. 2010b. A fluid limit for an overloaded X call center 

model via an averaging principle, working paper, Columbia University, NY. Available 

at: http: / /www. columbia.edu/~ww2040/allpapers. html 
[25] Salminen, P. and NORROS, I. 2001. On busy periods of the unbounded Brownian 

storage. Queueing Systems 39 317-333. 
[26] Rootzen, H. 1977. On the functional central limit theorem for Martingales. Zeit. 

Wahrscheinlichkeitsth. werv. Gebiete 38 199-210. 
[27] Talreja, R. and Whitt, W. 2009. Heavy-traffic limits for waiting times in many- 
server queues with abandonment. Ann. Appl. Prob. 19 (6) 2137-2175. 
[28] Whitt, W. 1980. Continuity of generalized semi-Markov processes. Math. Oper. Res. 

5 (4), 494-501. 

[29] Whitt, W. (1981). Comparing counting processes and queues. Adv. Appl. Prob. 13 
(1) 207-220. 

[30] Whitt, W. 1982. On the heavy-traffic limit theorem for GI/G/infinity queues. Ad- 
vances in Applied Probability 14 (1) 171-190. 

[31] Whitt, W. 1992. Asymptotic formulas for Markov processes with applications to 
simulation. Oper. Res. 40 (2) 279-291. 

[32] Whitt, W. 2002. Stochastic-Process Limits, New York: Springer. 

[33] Whitt, W. 2004. Efficiency-driven heavy-traffic approximations for many-server 
queues with abandonments. Management Sci. 50 (10), 1449-1461. 



Tech C110, 2145 Sheridan Rd. Department of Industrial Engineering and Operations Research, 

Evanston, IL 60208 Columbia University New York, New York 10027-6699, 

E-MAIL: ohad.perry@northwestcrn.edu E-MAIL: ww2040@columbia.edu 



